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Chapter 1 

Introduction 



One of the most recurrent philosophical questions in human thought has surely been: 
" What are things made of?\ 

People tried to answer this question during the times and found several and some- 
times eccentric solutions to the problem: from Greek philosophers to modern scien- 
tists many centuries of improvements and progress passed and the original question 
changed horizon time after time. 

Since the introduction of the galilean method of inquiry in science the search for the 
laws governing the world of the infinitely small building blocks of nature has been a 
great source of discoveries. The attempt to answer our initial question has increased 
dramatically the knowledge of mankind and has brought incredible applications in 
many sectors, up to be applied in everyday life of ordinary people (even if many 
people do not realize this, quantum mechanics has brilliant applications in many 
devices we use every day) . This kind of physics can look abstract and far away from 
ordinary needs, so that in non scientific environments the question about the sense 
of the construction of very expensive accelerators makes sometimes its appearance. 
However if we look at the past we can recognize that very important inventions were 
introduced on the basis of physical laws or phenomena that physicists discovered 
without any practical aim (no one worked out Quantum Mechanics to understand 
how to build a transistor nor anybody studied radiations, at the beginning, to radio- 
graph someone else's broken leg. . . ). The discoveries about the microscopic world 
that we are inquiring today will likely have a practical application in the future. 
Anyhow every new discovery about the world around us, at microscopic or macro- 
scopic or cosmic level, increases our knowledge of the universe and the nature and 
this constitutes a progress for our consciousness about what are things we can or 
cannot see around us, from the largest scales to the smallest distances. 
What is considered "small" is obviously a function of time: in the 19th century 
scientists began to introduce the concept of atom, an indivisible particle, building 
block of every state of matter. Just a few decades later, at the beginning of the 20th 
century scientists as Rutherford and Thompson showed that the indivisible atom was 
not so indivisible, but had its own inner components, nuclei and electrons. Moreover 
this microscopic world was described by new and unexpected laws, the Quantum 
Theory, sometimes in deep contradiction with our usual way of thinking. 
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In the 30s even nuclei began to show an inner structure, a bound state of neutrons 
and protons: little by little a new branch of physics was born, the Elementary Par- 
ticle Physics. To study the new features of elementary particles larger and larger 
energies needed to be reached and this was made possible by bigger and bigger 
instruments, accelerators and more and more advanced and refined detectors were 
needed too. At the beginning of this adventure, a powerful accelerator could be 
safely laid on a laboratory desk, nowadays people build accelerators of several kilo- 
meters of diameter wide to reach the astonishingly high energies necessary to go 
further in the exploration of small distances. 

As time went by the number and the content of discoveries concerning the world of 
extremely small distances began to separate from everyday life: before the discovery 
of the muons this branch of physics was studying particles involved in the matter we 
can touch every day. After that High Energy Physics began to study a world which 
exist over scales of time of the order of few microseconds at first, then nanoseconds 
and now even smaller times. This is a world with does not exist in ordinary life. 
A larger and larger number of short-living particles appear: this was a great puzzle 
for physicists trying to understand why nature was composed by so many different 
building blocks. Even if sometimes they do not admit, theoretical physicists are 
attracted by the idea that an ultimate theory of nature (both by elementary con- 
stituents and fundamental interactions) should be beautiful, though this assumption 
is surely difficult to define within a scientific framework. We use to think that beau- 
tiful in science is something simple, symmetric, which needs the smallest number 
of assumptions and ingredients. And a large number of building blocks is not a 
beautiful feature for a theory! 

The modern history of High Energy Physics begins in the 50s of the last century 
with Hofstadter's experiments at Stanford, which demonstrated that protons are 
extended objects and measured their form factors. In the following years a linear 
collider was built at Stanford and the experiments performed at SLAC (Stanford 
Linear Accelerator Center) brought to the resolution of proton components and the 
introduction of quarks. The hypothesis that hadrons are composed by point-like 
building blocks was introduced by Gell-Mann and Zweig, from spectroscopical ob- 
servations, and by Feynman and Bjorken, to explain the so called Bjorken scaling. In 
the same years the Standard Model was codified by Glashow, Weinberg and Salam. 
The latter is surely one of the most advanced goal reached by scientists in the 20th 
century: it was able to describe unexplained effects, to predict new phenomena. In 
practice all the quantitative predictions in High Energy Physics made by the Stan- 
dard Model are correct within the limit of experimental errors. 
Is this the final answer? The existence of three families of particles (even if only the 
first one is involved in ordinary matter) and the known fundamental forces? Many 
problems are still open as it is discussed in the following chapter. 
Physicists are not fully satisfied with this answer: they would like to include gravity 
in this picture, to understand where particle masses come from, why they interact 
in that way, why the families are just three,. . . . 

This is the reason why the main problem of this kind of physics nowadays is how 
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to go beyond the Standard Model: however the still correct predictions of the the- 
ory complicates the game, there are not data in sensational disagreement with the 
theoretical prediction, even if the discovery of a mass for neutrinos could open some 
door. 

In the last years, due to the progress made in understanding the strong and weak 
interactions, a new topic arose in High Energy Physics: the high precision study of 
the properties of the beauty quark, b physics. This quark was discovered in 1977 
and its role inside the Standard Model is peculiar: for contingent phenomenological 
reasons, in b physics many parameters very important to test the Standard Model 
with a high level of precision are involved. Moreover b physics can be very sensitive 
to effects due to phenomena that cannot be described by the Standard Model, the 
so called New Physics, a theory beyond the Standard Model which has not yet a 
precise form. 

Finally in the decays of the b quark the perturbation theory for the strong interac- 
tions can be applied, as we have done, since the property of asymptotic freedom, 
discovered by Wilczek, Politzer and Gross in 1974, implies that 2£ < 1, while for 
the decays of lightest quarks in general this is not true. 

A large program of measurements has begun in the past years, the final goal being 
a determination of the parameters concerning the physics of the b quark with high 
precision, to constraint the Standard Model and conclude if its predictions are fully 
compatible with experiments. 

This program is mainly based on the construction of the so called b factories, Belle 
and BaBar, accelerators built specifically to study the properties of the b quark. 
Also CLEO gave important results concerning this topic. In these years they were 
dedicated to the measurements of decay rates of heavy hadrons, lifetimes, branching 
ratios, parameters of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. In particu- 
lar the experiments gave the most precise measure of the elements involved in the 
unitary triangle and the first experimental observation of CP violation in b physics 
(with a measurement of the associated phase of the CKM matrix). 
At the same time new theoretical tools, based on the Standard Model, were intro- 
duced to face in the proper way the physics of this quark: in particular an effective 
theory, called Heavy Quark Effective Theory (HQET), found a large application in 
theoretical predictions, being able both to simplify the full theory and to reproduce 
its dynamics. 

This is the horizon where we decided to move our studies and inquiries: the specifi- 
cal motivations about the calculations and the analysis performed are in chapter 
here let us just introduce the topic of this thesis. 

In the framework of the b physics we decided to study, from a theoretical point of 
view, a rare process of decay of this quark, the transition b — > 57. This process has 
been widely studied in the past years because physicists hoped to see clear signals 
of new physics: this hope was frustrated by the good agreement of the theoreti- 
cal predictions with experimental data (even if for a few times a disagreement was 
observed). We focused our calculations on the transverse momentum distribution 
of the strange quark with respect to the direction defined by photon flight. The 
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kinematics of the process is introduced in chapter 

The main topic of this thesis is an application of the technique of resummation to 
this particular channel and an accurate and complete evaluation of the strong cor- 
rections. In particular a full 0(as) calculation and the structure of large logarithms 
are the crucial ingredients to reach this goal. 

The resummation of large logarithms is a technique introduced in QCD about 25 
years ago, to improve the perturbative expansion in regions of the phase space for 
processes where strong interactions are involved, for example strong radiative cor- 
rections to e + e~ — ► 2 jets or in semi-inclusive distributions in Drell-Yan. 
This technique has been developed and applied to processes where only light quarks 
were involved, since at the time accelerators were not dedicated to the intensive 
study of heavy quarks (they were just discovered). Here we will apply this tech- 
nique to a process were an heavy quark is present and this will change some of the 
dynamical features. 

Another important point we will argue is the reliability of perturbative QCD at 
relatively low energies: in fact perturbative QCD has had brilliant confirmations in 
high energy processes at the scale of the intermediate boson masses, while in this 
case the energy scale is a factor 20 smaller (Mz ~ 20mfe) and the strong coupling 
constant is about twice (as(Mz) ~ 0.12 ~ 2as{mb)). 

These problems will be discussed during the text, where the explicit calculations are 
shown. 

This thesis is organized as follows. 

In Chapter El aspects of the Standard Model are briefly described and in particular 
the most important features of QCD, involved in our following discussions, are out- 
lined. 

In Chapter El a brief review about the specifical process b — > is presented, above 
all a description of the effective hamiltonian used to disentangle the dynamics of 
this transition is given. 

Chapter H] contains the motivations for our work, the reasons why we decided to 
apply the technique of resummation to this particular process and against what 
background is set our analysis. In this chapter the kinematics of the process is also 
outlined. 

From Chapter El on the explicit calculations are shown: there the large logarithms 
appearing near the limit of the phase space are resummed in the impact parameter 
space, according to the method of resumming transverse momentum in an auxiliary 
space. 

In following Chapter El we take into account the singularities of the resummed distri- 
bution and we compare them with the analogous ones of another distribution con- 
cerning this process, the threshold distribution. This provides informations about 
non perturbative effects. 

Chapter [7| completes the calculation, with the evaluation of constants or regular 
terms which are not resummed and which do not show a logarithmic enhancement. 
Chapter El contains preliminary results about the effects of the introduction of a 
mass in the final state, for example for the strange quark: this is a necessary step 
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to deal with with another important transition of the b quark, namely b — > c. 
Finally Chapter |H1 a brief summary of the contents of thesis is presented with an 
outlook of future improvements. 

Appendices, which deepen some topic met during the text, are given at the end. 

Some of the results obtained during the work for my thesis have been already pre- 
sented in Pig. 
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Part I 

Physics of the b Quark and the 
Standard Model 
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Chapter 2 

Standard Model 



2.1 The Standard Model of Elementary Particles 

The Standard Model is nowadays the framework describing interactions among ele- 
mentary particles in processes occurring in a range of energies up to about 1 TeV. 
It is a non abelian and renormalizable gauge theory based on the symmetry group 

Q SM = SU(3) C ® SU(2) L ® U{l) Y . (2.1) 

The sector corresponding to the group SU(3)c describes strong interactions through 
Quantum Chromodynamics (QCD), which will be discussed in this chapter from sec- 
tion (|2.Hj) on, while the sector corresponding to the group SU{2)l®U{1)y describes 
the unification of weak and electromagnetic interactions through the Glashow— 
Weinberg-Salam model [SlHj, briefly described in section (|2.2p . 
The original symmetry (|2.1|) is broken by a non vanishing vacuum expectation value 
in the theory for a scalar field, which gives mass to the gauge bosons of the elec- 
troweak sector and to the fermions. 

The main details of theory are discussed in the next section, here we will just make 
some observations about the model: 

• up to now experiments performed in the laboratories all over the world confirm 
the predictions of the Standard Model: this statement could seem too optimistic 
and strong, but, as a matter of fact, there are not observations or calculations 
which are within the errors violated j^]; 

• the recent results obtained in agreement with hypothesis of neutrinos oscilla- 
tions seem to be a highly non trivial check of the Standard Model: we are not 
currently able to state wether this phenomenon can be brought back to the 
Standard Model or is a signal of new physics, leading to a theory beyond. In 
principle a mass for the neutrinos could be introduced in the Standard Model, 
as discussed later. Other searches of new physics (studying rare decays of the 
beauty quark, measuring CKM matrix elements and so on) have not given any 
result. Up to now no signals of supersymmetry, extra dimensions, technicolor 
have been found; 
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• the Higgs boson has not yet been detected: one of the most important goals 
of future accelerators is its discovery. If the Higgs boson does exist it will be 
possibly detected at LHC: the most disappointing eventuality would be the 
discovery of a single Higgs particle as predicted by the Standard Model, since 
this will not give any information about extensions of the Standard Model; 

• though its success the Standard Model is not believed to be a fundamental 
theory for several (and good) reasons: the most relevant one is that it does 
not describe gravitational interactions. At present a quantum description of 
gravity is lacking; 

• moreover it depends on many free parameters (34) which are fixed by experi- 
ments and cannot be calculated from the theory: a final theory should be able 
to explain why the mass of particles have their values, why coupling behave 
that way and so on and the answer to this question is still lacking; 

• the interactions described in the Standard Model are not really unified, because 
they depends on different and not related couplings. Maybe they are unified 
at energies much larger than the ones reachable by present machines. How- 
ever without a supersymmetric model the couplings of the interactions in the 
Standard Model do not unify; 

• furthermore the Standard Model does not explain correctly the asymmetry 
between matter and anti-matter. Inside the theory the CP violation (in the 
last years confirmed by Belle and BaBar) could introduce such an asymmetry, 
but the quantitative prevision is largely incorrect; 

• other issues (such as the hierarchy problem) suggest that the theory should 
be valid up to a new physics scale A, where supersymmetry should make its 
appearance. This is the hope of many high energy physicists: in my opin- 
ion, lacking strong experimental evidences that violate the theory, it is very 
complicated to put forward the correct guess. 

By the light of these observations we can conclude that the Standard Model is 
an excellent parametrization of particle physics and has been able to explain high 
energy processes studied in experiments performed up to now. 
The theoretical issues, as the ones mentioned above, suggest that this is not the 
final theory of particle physics and fundamental interactions: a theory beyond the 
Standard Model is the Holy Graal of modern particle physicists, but its elaboration 
collide with the lack of any signal of new physics. 

In this chapter the main features of the Standard Model will be recalled, with a 
particular emphasis to its strong sector. 
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2.2 Glashow-Weinberg-Salam Model for Electroweak Inter- 
actions 

In this section the main properties of the Glashow-Weinberg-Salam model [Sill], the 
theory of electroweak interactions, are briefly recalled. 

This theory was put forward in the 60s, as a description of the unification between 
weak and electromagnetic interactions: several attempts were tried to introduce a 
theory of massive bosons for weak interactions, but this kind of models turn out to 
be non renormalizable. Moreover physicists realized that also the photon had to be 
included in this description, to make the theory predictive 1 . 

Experimental observations suggested that weak interactions does not conserve parity 
jHj and in particular that only left handed fermionic fields interacts weakly: this 
drove to build a lagrangian with massless fermionic fields, with independent left- 
handed and right-handed components: 

1 — 7c; 

*r = (2-2) 

The model is based on the hypothesis that fermionic left-handed fields are arranged 
in isospin doublets, while fermionic right-handed fields are singlets. 
Fermionic elementary fields, leptons and quarks, are rearranged into three families, 
each containing a lepton doublet and a quark doublet of left-handed fields: 

^\ ( \ / V T 
/I J \ T 

(2-3) 

u \ I c \ ft 



d J I s j \ b 

Right-handed fields are singlet of isospin: 

£r, Hr, ur, d R , ... (2.4) 

In order to consider also electromagnetic interactions in the model a new quantum 
number, the hypercharge Y, was introduced. The hypercharge of a particle is defined 
as twice the difference of its electric charge and the third component of its isospin: 

Y = 2(Q-T 3 ). (2.5) 

The gauge group of symmetry of the model is: 

yaws = SU{2) L ®U{1) Y . (2.6) 

Such a group has 4 generators, corresponding to the gauge bosons of the theory, 
A a (x), with a = 1, 2, 3, and B^(x): in order to get the invariance under the action 
of the group, a covariant derivative is introduced as 

D» = ffl * - igT a ■ A" - ig'B". (2.7) 

lr The renormalizability of the theory was shown by 't-Hooft and Vcltman in |7| 
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The kinematical term for the gauge bosons is introduced through the strength tensor: 



A,-,, = -\g; u g^ - \f» v f»\ (2.8) 



defined as 



Ffx V = d^-dvB^. (2.9) 

The kinematical part of the lagrangian contains no mass terms for the gauge bosons: 
the experimental evidence that weak interactions has a small range of action , about 
10 _17 m, suggest that these interactions are mediated by a massive boson, whose mass 
is about 10 2 GeV. As it was demonstrated, a mass term introduced by hand would 
lead to a non renormalizable theory. The mass is given to the gauge bosons (and 
to the fermions) by the mechanism of spontaneous symmetry breaking, through the 
introduction of a scalar field $ which acquire a non zero vacuum expectation value, 
< $ > ^ 0. 

One introduces the scalar doublet jBJ 



*v> = \\ )• < 210) 

and in the lagrangian the term 

Cmggs = \D^? + - A(0*0) 2 . (2.11) 

In this way the non vanishing vacuum expectation value reads 

and, by using the gauge symmetry of the theory it can be rearranged as 

^ ) = 7f(, + V))' (2 - 13) 

In this way the scalar field can be expanded around its minimum < <\> >= v/y/2, 
introducing the scalar field 77 (x) which represents the Higgs boson. 
The covariant derivative in Cm ggs produces terms of interactions between the gauge 
bosons and the scalar doublet: once the symmetry is broken, by giving a non van- 
ishing vacuum expectation value to 0, mass terms for the gauge bosons arise: 

j(gA$ - g'B»Y + V ^g\A x + %Atf(A x - iA 2 ),. (2.14) 

beside terms of interactions between the gauge bosons and the Higgs boson, rj(x). 
The terms in ()2.14|) can be diagonalized in order to interpret them as mass terms in 



2.2. GLASHOW-WEINBERG-SALAM MODEL FOR ELECTROWEAK INTERACTIONS 19 



the lagrangian for the gauge bosons: this can be accomplished defining 

W ± = -=(A 1 ±iA 2 ) 
v2 

Z = — B sin 9 W + A 3 cos 6*^ 

7 = B cos 6 W + A 3 sm6 w , (2.15) 

sin% = ^=L=. (2.16) 
V9 + 9 

This implies that bosons W ± and Z° acquire a mass, while the field 7, which can 
be interpreted as the photon, remains massless: 



where 



gv y/g 2 + ff' 2 
m 7 = 0, m w = m z = . (2.17) 

The electric charge can be identified with the following combination of parameters 
of the theory 

e = gsin6 w . (2.18) 



Moreover, from a comparison with the Fermi theory of weak interactions, one can 
state that 

'F 



g 2 Gi 



2 



V2~ 



[2.19) 



At the end the original symmetry is broken into 

SU(2) L ®U(l) Y ^U(l) em . (2.20) 

The interaction of the fermionic fields (leptons and quarks) is introduced as usual in 
the construction of a gauge theory, considering the Dirac lagrangian and imposing 
the invariance under the symmetry group: this implies that the ordinary derivative 
have to be replaced by the covariant derivative ()2.7|) . 

The mass of fermions cannot be introduced by hand, because the massive term in 
the Dirac lagrangian would violate the symmetry SU(2)l, so that 

C ferm =^ L (i^D^ L + ^ R {i^D^ R . (2.21) 

Fermions acquire a mass through the symmetry breaking, by introducing in the la- 
grangian a Yukawa type interactions between (left-handed and right-handed) fermion 
fields and the scalar field, that is terms such as 

^Yukawa = 9f{Le R (f) + e R L(f)), (2.22) 

After symmetry breaking a mass for fermions arise: 

ZW = E^7/, m/ = M. (2.23) 
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where the summation runs over the different fermions (possibly except neutrinos). 
Moreover terms of interaction among the Higgs boson and the fermions 



showing that the Higgs boson couples to fermions with a factor proportional to their 
mass. 

Unfortunately the parameters gf cannot be calculated from the theory, this fact also 
indicating that the Standard Model cannot be considered a complete description of 
the particle world. 

A different discussion for neutrino masses can be made: one should introduce right- 
handed neutrinos, but being the neutrino a chargeless particle it interacts neither 
through electromagnetic interactions, and the right-handed part nor through weak 
interactions, so that the introduction of a right-handed neutrino seems not to be a 
necessary step in the Standard Model. Recent measurements indicates that neutri- 
nos do have a mass: this implies that or a right-handed neutrino is introduced in 
the theory, as done for other fermions, or other mechanisms of generation of a mass 
are necessary, which would eventually point to new physics effects. At present the 
problem of neutrino masses and its implications in the search of new physics is still 
open. 

The last point to deal with to get the lagrangian of the Standard Model is the mixing 
between quarks. An analogous phenomenon for charged leptons can be introduced, 
but it has been poorly studied and is at the moment less relevant, while neutrino 
mixing is one of the major topics, related to the measurement of their masses through 
oscillations. Writing down the lagrangian of the Standard Model one realizes that 
quark fields interacting with gauge bosons and quark fields mass eigenvectors are 
not forced to be the same, because they appear in different contexts. This allows 
to introduce a unitary matrix Vqkm, the Cabibbo-Kobayashi-Maskawa matrix 2 [5], 
which describes the relation between the eigenvectors of the interaction hamiltonian 
q' and the mass eigenvectors q: 



The choice of down type quarks for the mixing is purely a matter of convention: 
a quark mixing using up type quarks could be considered and physical prediction 
would not be affected by this choice. 

The existence of this matrix explains transitions between different families, which 
would not be allowed without it: CKM matrix elements are in practice a measure 
of the strength of interactions between two quarks of different flavour. 




.f 



(2.24) 




(2.25) 



2 Cabibbo at first introduced a mixing between the two first generations of quarks. Kobayashi and Maskawa 
extended it to the heaviest one. 
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The CKM matrix can be parametrized and written in different ways, for example: 

( V ud V us V uh \ ( 1-f A AA 3 (p-^)\ 

Vckm = V cd V cs V ch - -A 1 - f AX 3 

\ V td V ts V tb ) \ A\\l -p-irj) -A)? 1 / 

(2.26) 

The last term in ()2.26j) is the Wolfenstein parametrization [Hlj: it is an approximated 
form, valid up to 0(A 4 ), where A is the Cabibbo angle, r] parametrizes the complex 
part of the matrix and it is a measure of CP violation in the Standard Model. In 
the last years, at the b-factories, a large effort has been made to measure, with 
high precision, CKM matrix elements, above all the smallest ones, in order to verify 
the unitarity of the matrix. A violation of unitarity could be a possible signal that 
something is incomplete in the theory. In particular most of the experiments tried 
to verify the relation 

v ud v: b + v cd v; b + v td v; b = o (2.27) 

also known as unitarity triangle 3 . Latest experimental results give the following 
values for the parameter in the Wolfenstein parametrization: 

A = 0.2279 ±0.0032 (2.28) 

A = 0.768- 0.824 (2.29) 

p = 0.118 - 0.273 (2.30) 

Tj = 0.305 - 0.393. (2.31) 

3 Unitarity triangles are in general three and corresponds to the constraint of unitarity for the elements of the 
matrix Vckm- However, for phenomenological reasons, the other two triangles are more difficult to test, so that 
the unitarity triangle for antonomasia is the one in 12.271 . 
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These values are consistent with the unitarity of the CKM matrix, as shown in the 
figure (j2.1|) [11. Once again the Standard Model seems to withstand experimental 
tests. 

2.3 Quantum Chromodynamics 

Hadrons are particles interacting through strong interactions: since the 30s-40s, from 
Yukawa's theory of strong forces on, they were studied to find a coherent theory of 
this kind of forces. Theoretical attempts have clashed with the large coupling con- 
stant associated to these interactions which would have precluded any perturbative 
expansion. 

This difficulty was overcome in the 70s with the discovery of Quantum Chromo- 
dynamics (QCD). QCD is a non abelian gauge theory introducing a new quantum 
number, called color, which exists in three different forms: red, blue and green. The 
fundamental fields of the theory are the quarks , which carry the charge of color 
and interact with each other by exchanging the carriers of the strong force, called 
gluons. 

Experiments suggest that colored particles cannot be observed as free states in na- 
ture and to explain this evidence an additional property, the confinement, has been 
postulated. 

In this picture, colored particles are confined into hadrons and they are combined 
in such a way that the color is screened outside the hadron and does not show up. 
A description of the confinement is not yet possible: in particular we cannot solve 
completely the theory, starting from the lagrangian of QCD, written in term of fun- 
damental fields (quarks and gluons), to explain the spectrum of the bound states 
observed in nature, i.e. the hadrons. 

In spite of the lacking of this complete solution, QCD can explain in great detail 
strong interactions in high energy experiments. 

Although strong interactions involve a large coupling constant at low energies, the 
application of renormalization group equations shows that, at higher energies, the 
coupling constant becomes rather small and a perturbative approach is viable. 
This occurs because 

a s {Q 2 ) < 1 for Q 2 > A 2 QCD (2.32) 

where Aq C d is the typical scale of strong interactions, Aq C d ~ 200 MeV and Q 
is the hard scale of the process, that is the highest energy scale where the process 
takes place. 

At energies comparable with Aqcd the coupling constant blows up and any pertur- 
bative approach is barred. This behaviour can qualitatively explain the origin of 
confinement: when quarks separate (large distances correspond to small exchanged 
momenta), the interaction between them becomes stronger and stronger, in such a 
way that they can never be further than a distance corresponding to the typical 
scale of QCD. 

During the times different techniques were developed in field theory to deal with the 
non perturbative region of strong interactions, such as for example QCD sum rules 
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and lattice QCD. 

On the contrary the perturbative part has been widely studied and tested in the last 
twenty five years and perturbative QCD has become a precision tool to deal with 
strong interactions in high energy physics. 

In the next sections some of the most important features of QCD will be briefly 
recalled, above all those properties which are needed further down. 



2.4 The Parton Model 



In the 60s experiments performed at the Stanford Linear Accelerator Center (SLAC), 
showing the so called Bjorken scaling |T2J , and spectroscopical observations suggest 
that mesons (bosonic hadrons) and baryons (fermionic hadrons) are composed by 
building blocks which, at high energies, interact with external currents as free par- 
ticles. 

These constituents were called partons, the name introduced by Feynman, or quarks, 
as they were called by Gell-Mann On theoretical basis also Zweig |T3] intro- 
duced a similar entity 4 . 

Experiments at SLAC, pointing out the partonic behaviour at high energies, con- 
cerned Deep Inelastic Scattering, that is the scattering of an electron over a proton 
with a large exchanged momentum (that is a momentum much larger than the 
proton mass). In terms of the Bjorken variable 

x = ~w~ q < 2 - 33 > 

where g M is the exchanged momentum and P M the proton momentum (usually taken 
in its rest frame, where P M = (M, 0)). In first approximation the structure functions 
f a {x), representing the probability to find a parton of type a in the hadron with a 
value x of the Bjorken variable, measured in the process do not depend on x (Bjorken 
scaling). As it was noticed, this is the typical situation occurring in point-like elastic 
scatterings: this behaviour can be explained assuming that the electron is scattered 
by a point-like constituent inside the proton, suggesting that at high energies hadrons 
can be thought as composed by partons, which behave as free particles. 
At the beginning of the 60s, to explain hadron spectroscopy, three types of quarks 
were introduced: up, down and strange. In the following years other three quarks 
were discovered: the charm quark in 1974, bottom or beauty in 1977 and , finally, 
the top quark in 1994. 

The first three quarks have masses negligible at high energies and for this reason 
they are usually referred as light quarks. The three quarks discovered more recently 
are referred as heavy quark, because their mass is in general relevant. 

m c ~ 1.5GeV, 
m b ~ 5GeV, 

m t ~ 175GeV. (2.34) 



*Zweig introduced the name aces. His paper was not published. 
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This is true for the top and, in many processes (and in particular the ones we are 
interested in), for the beauty. The mass of the charm quark is relevant in many pro- 
cesses, but the application of perturbative QCD at charm mass energies is doubtful, 
since the coupling constant is quite large as{m c ) = g 2 /4:TT ~ 0.35. 
The parton model was introduced by Feynman and Bjorken for light quarks, assum- 
ing that each parton in the hadron carries a fraction of the total quadrimomentum of 
the hadron itself. Defining f a (z) the probability to find a parton of type a, carrying 
a momentum fraction z, the cross section da for the Deep Inelastic Scattering can 
be written as 

da = dz f a (z) da a (z) (2.35) 

a J 

where da is the partonic cross section, that is the cross section for the scattering of 
an electron over a quark a. 

At first the parton model was a naive picture which did not take into account strong 
corrections: an improved parton model can be consistently introduced considering 
QCD corrections, which explain for example Bjorken scaling violations. 
Eq. (l2.35j) shows a paradigmatic situation in QCD, the interplay between non pertur- 
bative and perturbative terms, which contribute together to the result, but can be 
obtained in different ways: da a (z) being a point-like object can be in principle com- 
puted in perturbation theory, while the parton distribution functions (pdf), being 
log distance objects, cannot be calculated and have to be extracted by experimental 
data. However the theory allows to calculate perturbatively the evolution of the 
pdf 's as a function of the scale. 



2.5 The Lagrangian of Quantum Chromodynamics 

QCD is a non abelian gauge theory, whose group of symmetry is SU(3) C . The index 
c indicates that the quantum number is the color. 

The Lagrangian of QCD is obtained according to Yang-Mills theories, by requiring 
a local invariance under the group SU(3) C 5 . The construction of the Lagrangian 
involves the definition of a covariant derivative 

{D») ah = d»8 ah + ig s {t c A^) ab (2.36) 

where A^(x) are gluons, the gauge fields of the theory, and t a are matrices of the 
fundamental representation of 57/(3), having the properties 

[t A ,t B ] =if ABC t c (2.37) 

being f ABC the structure constants of SU(3). 

The most common choice is provided by the eight Gell-Mann matrices, hermitean 
and traceless, normalized as 

Tr t A t B = l -b AB (2.38) 

5 Since now the suffix c will be neglected. 
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Important relations satisfied by the colour matrices are 

Y< t abti = C F 5 ac with C F = \. (2.39) 

A 6 

For the adjoint representation it holds 

Tr T A T B = C A 6 AB with C A = 3 (2.40) 

The lagrangian density describing the interaction between fermionic fields q(x) 
(quark) and gauge bosons (gluons), locally invariant under SU(3) turns out to be: 

Cqcd = -\FZ v F a » v + U*P-rn) ab q b (2.41) 

flavours 

where F^ u is the strength tensor defined as 

F; v = [d,A a v - d u A« - g s f abc AlAl\ . (2.42) 

The indices a, b, c run over the eight degrees of freedom of the gluon field. 
Feynman rules deriving from this lagrangian are contained in appendix (jXjl . 



2.6 Evolution of the Coupling Constant and Asymptotic 
Freedom 



Let us consider an observable 0(Q 2 ), where Q is an energy scale much larger than 
every other energy parameter involved in the process, for example the masses of the 
quarks, Q 2 ^> m 2 . In general a prediction about the value of this observable can be 
performed in perturbation theory, that is as an expansion in the coupling constant 



a. 



9± 
4tt' 



(2.43) 



provided that as is small enough to justify this approach. Obviously let us consider 
an observable where only strong corrections need to be calculated. 
In order to remove ultraviolet divergences, arising to every order of the perturbative 
expansion, the theory of renormalization has to be applied to the observable 0(Q 2 ), 
by introducing a substraction point \i. Now, having introduced a second energy 
scale, the observable O will depend on the ratio Q 2 / '/i 2 . 

However the renormalization scale \x is arbitrary and, according to the theory of 
renormalization, the observable (9(Q 2 //i 2 ) has to satisfy the Callan-Symanzik equa- 
tion 6 (renormalization group equation): 



- d n(® 2 \ 



dfi 2 



_d_ 



+ /r 



,das d 
d[i 2 da. 



s 



O 2 



0. 



(2.44) 



3 'We are neglecting the masses and this will lead to simpler equations. 
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Defining the (3 function as 

= /i 2 || (2.45) 
and defining t = log(Q 2 //i 2 ), the equation (|2.44J) can be written as 

O(e t ,a s ) = 0. (2.46) 



+ /3(a s )- C 



dt das 
An implicit solution of this equation is given by 

0(n 2 = Q\ a s {Q 2 )) = 0(1, a s {Q 2 )\ (2.47) 

where the new function as{Q 2 ), the running coupling constant , is defined according 
to the property 

" asiQ2 ^ da , (Q 2 



Pi = 10 2 S (2-50) 



, log Pj • ( 2 - 48 ) 

Let us remark that in this way the whole dependence on the scale is absorbed in 
the running coupling constant. 

The function f3(as) admits a perturbative expansion in the form 

P(a s ) = -Poa 2 s - - f3 2 a% - . . . (2.49) 

and substituting ()2.49|) into ()2.45j) . the running of the coupling can be calculated in 
perturbation theory. 

The coefficients (3 n can be calculated from higher order corrections to the bare 
vertices of the theory: at present they are known up to fourth order. For our 
purposes only /? and (3\ will be relevant: 

33 - 2n f 
153 - 19n f 

where n/ is the number of the active flavours at the energies where the process takes 
place. 

The leading order solution of the equation ()2.45|) reads 

In this way we can extract the value of the coupling constant at a scale Q, known its 
value at Qo. The most common choice is to consider as a reference as(M|o) which 
is a parameter of the Standard Model known with very high precision. 

a 5 (M| ) = 0.11720 ± 0.00017 (2.52) 

Another choice to express the solution of ()2.45|) is possible: one can introduce a 
parameter Aqcd m such a way that ()2.51|) takes the form 

Air 

as ^ = 7H (H2,K2 V ( 2 - 53 ) 
P \og(Q 2 /A 2 QCD ) 
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In this way the running coupling constant is expressed in terms of one single param- 
eter and it is evident that it blows up when the energy scale is set equal to Aqcd- 
The exact value of Aqcd strongly depends on its precise definition, but it can be 
considered of the order of Aqcd ~ 200 MeV. Roughly speaking, for energies of the 
order of a few hundred MeV, the perturbative approach for strong interactions is 
not reliable because as becomes too large. Let us notice that this is the typical scale 
of masses of light hadrons and this suggests that the large growth of the running 
coupling at low energies is an ingredient necessary to explain phenomena such as 
the confinement of the quarks inside the hadrons. 

On the contrary when the scale of the process becomes very large, the running 
coupling constant vanishes: 



This property is known as asymptotic freedom and was discovered by Wilczek, 
Politzer and Gross in 1974 ^2]- Asymptotic freedom is a fundamental property 



• it justifies the perturbative approach for energies much larger than Aqcd'- since 
as becomes small enough, one can expand any observable in powers of this 
parameter; 

• it explains why the parton model is a successful picture of the structure of an 
hadron: for large Q 2 the partons in the hadrons behave as non interacting par- 
ticles in first approximation, because the coupling constant is small. Through 
radiative corrections one can calculate violations to this behaviour. 

Let us observe that the existence of asymptotic freedom is related to the sign of (3$: 
as long as the number of flavour is n/ < 33/2 (as it happens in nature), the first 
coefficient /3q is positive, that is the first term in the expansion ()2.49j) is negative. 
A step-by-step solution of equation ()2.44j) would show that as{Q 2 ) becomes smaller 
and smaller as the energy becomes larger because of the sign of (3q. 
This behaviour is opposite with respect to QED, where (3q = — (according to the 
conventions of eq. ()2.49|) ) [ Tfij . 

The physical origin of this difference arises from the non abelian nature of QCD: 
being a non abelian theory gauge bosons carry the quantum number of the symmetry 
group, the color, contrarily to what follows in QED, where photons do not carry 
any electric charge. This implies that gluons can interact with each other, as shown 
by Feynman rules in appendix El. In QED the coupling constant decreases at large 
distances (small energies) and this effect is naively explained as a consequence of 
vacuum polarization by electron-positron pairs, which screen the electric charge. 
In QCD this effect related to quark interactions is overwhelmed by an anti-screen 
effect due to gluons self-interactions, which causes a growth of the coupling constant 
at large distances. 

Just to conclude this section let us notice that in next applications, instead of ()2.53|) . 




oo. 




of QCD: 
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we will use the following expression for as{Q 2 )'- 

(3, log log {Q 2 /A 2 QCD ) 
Po log(QVA^) 



as(Q) 



A, log (QVA 2 QCD ) 



(2.55) 



which represents the running coupling constant calculated with next-to-leading ac- 
curacy. 

2.7 Infrared Divergencies 

Infrared divergencies and their cancellation, which is a crucial test of the consistency 
of a massless theory as QCD and QED, will be one of the main topics of this thesis. 
Let us consider a generic process involving at least a quark as an external leg and 
assume that the quark mass can be neglected with respect to the other relevant 
energy scales. 

Strong corrections to this process involve gluon bremsstrahlung: for example, for 
one real gluon emission, a straightforward calculation shows that the emission rate 
turns out to be proportional to the term 

r 1 du r 1 de 

r oc / — / - — , 2.56 

Jo v Jo (1-COS0) 

where u is an adimensional and unitary variable, defined as 



E gluon , ... , 

W = ~pMAX ( 2 - 57 ) 

gluon 



and 9 is the angle between the gluon and the quark. 

The expression in (|2.5fi|) is manifestly divergent and, in particular, it shows two 
types of singularities: 

• when the gluon momentum vanishes 

u -> (2.58) 

the so called soft singularity arises. It originates from the massless nature of 
the gluon and can be regularized by giving to it a fictious mass A; 

• when the angle of emission vanishes 

9^0 (2.59) 

the singularity is called collinear. It can be regularized by giving a mass m to 
the quark. 

Soft and collinear singularities will be referred with the comprehensive term of in- 
frared divergencies (or singularities). 

In general infrared singularities arise in theories when a massless field is present 
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(soft divergence), the gluon in QCD, or when it couples to another massless field 
(collinear divergence), for example to a quark [TTj . 
With the regularization indicated above, (|2.5fij) reads 7 : 



After two simple integrations one gets that the singularities are parametrized in 
terms such as log A and logm. 

This regularization is not the one we will use in explicit calculation: we will use 
dimensional regularization ^H] , which is more elegant and theoretically useful, even 
if it hides the origin of the singularities discussed above. In practice one calculates 
integrals in N = 4 + e dimensions, instead of 4 dimensions: in this number of dimen- 
sions infrared divergent integrals becomes regular and singularities are parametrized 
as poles in the regulator e. Once divergencies are removed in such a way, one can 
consider the limit e — > 0, coming back to a physical number of dimensions. 
In this regularization scheme the equation ()2.56|) takes the form 



More details about this topic will be discussed in chapter where explicit calculation 
will be shown. 

2.8 Cancellation of Infrared Divergencies 

Obviously meaningful physical results cannot be affected by singularities: somehow 
infrared divergencies have to be cancelled. 

At first let us consider the case of inclusive quantities, such as for example the total 
rate of a scattering or the total width of a decay. 

The non abelian nature of QCD complicates the question: let us recall at first the 
solution for QED, which is an abelian theory. 

In QED it is known from a very long time that infrared divergencies arise in photon 
brehmsstrahlung processes and in radiative corrections in general. Their cancella- 
tion is assured by Bloch-Nordsieck theorem ^Hj: it states that infrared singularities 
cancel if we sum over all the degenerate final states. In practice, order by order in 
perturbation theory, one has to sum all the real and virtual diagrams describing the 
same process: for example, for one photon emission one has to take into account 
photon brehmasstrahlung and virtual one loop diagrams with no photon emission. 
From a physical point of view, the explanation of this theorem is rather intuitive: for 
soft emission , for example, one should consider that every real detector has a finite 
resolution and cannot distinguish between an electron and an electron plus a soft 
photon 8 . In practice the bare state of electron does not exist: a physical electron is 
always surrounded by a cloud of soft photons. 

7 For simplicity let m be an adimensional variable proportional to the quark mass. 
8 In practice a photon with an energy lower than the energy resolution of the detector. 






(2.61) 
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Analogous considerations can be made for collinear singularities: every detector can- 
not resolve particles in a arbitrarily small angular cone, so that a summation over 
the state of electron and electron plus a cloud of collinear photon is required. 
However in non abelian theories the Bloch-Nordsieck theorem is in general violated, 
as a direct consequence of the non abelian nature of the theory fH\ I24j . 
The cancellation of infrared singularities in non abelian theories is assured, on the ba- 
sis of general arguments, by the Kinoshita-Lee-Nauenberg theorem, which provides 
a generalization of Bloch-Nordsieck assertions. The main results of this important 
theorem are recalled in the next section. 

2.8.1 Kinoshita-Lee-Nauenberg Theorem 

According to the Kinoshita-Lee-Nauenberg theorem [221 E3| the cancellation is at- 
tained when a summation over final degenerate states (like in Bloch-Nordsieck the- 
orem) and initial degenerate states (due to the non abelian nature of the theory) is 
performed. In particular, in Kinoshita's work [22], the relation between the Feyn- 



man diagrams involved in the process and the cancellation of mass (another usual 
name for collinear divergencies) and soft singularities is widely discussed. 
Let us recall the main results of |22]: let us consider a process at a given order 
of perturbation theory; let us call T the corresponding Feynman amplitude (see 
Fig. ([2.2Jl ). the total transition probability is proportional to J^^T-Tj, where the 

sum over indices is performed considering diagrams with the same final states. T/ is 
represented by a diagram obtained by Tj by reversing time, that is by the exchange 
of initial and final states. 

In this way T\Tj may be represented joining the final states of Tj and the initial 
states of T- (that is the final states of Tj): in order to distinguish it from a Feynman 
diagram let us draw a line intersecting the final states of Tj and Tj. 
Let us call T d the Feynman diagram obtained removing the cutting line and consider 
all the T\Tj which reduces to T d when the line is removed: by the optical theorem 




Figure 2.2: Example of two diagrams having the same final states. 



2.8. CANCELLATION OF INFRARED DIVERGENCIES 



31 



^2ijT^Tj, with the sum over the diagrams which reduce to T d when the cutting 
line is removed, is the absorbitive part A d of T d and it's called cut diagram (see 
Fig. (j2.3|) ). The total transition probability is the sum of cut diagrams involved in 




withoul; cuts rpd 



Ti + T 2 = l m r d 



Figure 2.3: Example of cut diagrams: Ti and T2 are two different cuts of the same diagram T d . 

the process. 

Let us now connect the initial states of T d with its final states (remember that they 
carry the same momentum), that is the initial states of Tj and Tj and let us represent 
this junction with a second cutting line. 

Removing the cutting lines we obtain a vacuum-to- vacuum transition T v : let us 
consider the set A of diagrams Aj which reduce to T v by removing the cutting lines 
(see Fig.flHJ)). 

A is called a double cut diagram and the total transition probability is a sum of As. 
Kinoshita's paper [22! shows that not only the total transition probability is free 




Ai A 2 A 3 



Figure 2.4: Example of double cut diagram: the amplitudes Ai,Aa and A3 belong to the same set A. 
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from infrared divergences but also every A doesn't show neither soft singularities 
nor collinear (or mass) ones. 

2.8.2 Cancellation of Infrared Singularities in QCD 

In principle, in QCD the summation over initial states would be required, according 
to what stated above. 

In many processes however Bloch-Nordsieck theorem can be applied, provided that 
a summation over initial and final color is performed; actually this is first of all a 
phenomenological requirement: because of confinement colored particles cannot be 
observed, so that a summation over all the possible color is natural, since only this 
superposition is observable. 

However some counterexamples have been found, where the only summation over 
colors is not sufficient: for examples in [2^ has been found that in Drell-Yan 



where X is every hadronic state, subleading soft divergencies arise at two loops and 
they are not cancelled by soft gluon emissions. If the degeneracy of the initial quark 
with the state of a quark plus a soft gluon is taken into account, this divergence 
is cancelled. This is nothing but an application of the Kinoshita-Lee-Nauenberg 
theorem. 

Let us underline that this is not necessary for our purposes and in our explicit 
calculation the simple summation over colors and final states will be enough. 

2.8.3 An Example: Incomplete Cancellation of Infrared Logarithms in 
the Electroweak Sector 

In the last few years, some authors j2"Ul |2"T] showed that Bloch-Nordsieck violations 
can be observed also in the electroweak sector of the Standard Model. This happens 
since the Glashow-Weinberg-Salam model is based on a non abelian gauge theory. 
It assumes particular characteristics because of some peculiarity of the theory. 
Apart the photon, whose behaviour in QED is discussed above and it is not inter- 
esting here, the gauge bosons of electroweak interactions, W ± and Z°, are massive. 
The mass provides a natural cut off: in this case infrared singularities are screened, 
but residual large logarithms, such as 



arise as in ()2.6()|) . though here the cut off mass is physical and not fictious. 
In |2"TH l2"T] authors noticed that these logarithmic terms do not cancel even in inclu- 
sive quantities: in QED and QCD large infrared logarithms appear in semi-inclusive 
distributions (they will be introduced in section (|2.10Jl ). but cancel in inclusive quan- 
tities, such as total rates. 

The appearance of residual logarithms in the Electroweak Model is a violation of 
the Bloch-Nordsieck theorem, due to the non abelian nature of the theory: in this 



qq -> 7X, 



(2.62) 




(2.63) 



2.9. EVOLUTION EQUATIONS 



33 



case a summation over initial weak charges is not required, contrarily to what per- 
formed in QCD with initial colors, because particles carrying weak isospin can exist 
as asymptotic states, for example as electrons or neutrinos. 
This implies that in processes such as 

e + e~ —>■ everything (2.64) 

terms like ()2.63|) remain. 

A cancellation of these terms would be achieved by summing the rate cr(e + e~ — > 
everything) to the rate cr(V e + e~ — > everything): this corresponds to a summation 
over initial state weak charges and, according to KLN theorem, provides the can- 
cellation of logarithms. Obviously there is not a compelling physical reason to do 
that. 

The presence of terms (|2.fiH|) does not spoil the theory: after all they are finite and 
not singular, even if they are a sort of shadow of the infrared singularity. 
However, at asymptotic energies, these terms can become large and enhance elec- 
troweak corrections, even if aw is quite small, making electroweak corrections com- 
parable to strong ones. If this effect does exist, it will be surely detectable at future 
Linear Colliders. 

2.9 Evolution Equations 

Evolutions equations are an important property of QCD dynamics |25| I26j. 
Let us consider multiple branching of partons from another one, for example multiple 
emissions of gluons from a quark or splitting of a gluon into a quark-antiquark pair 
and so on. As shown in picture ([2.5)1 . the parent parton, participating to an hard 




Figure 2.5: Evolution of a light quark outgoing from an hard process at the energy Q 2 . Through 
multiple gluon emissions it evolves from the virtuality t n to the virtuality to. 

process at a scale Q 2 , evolves emitting other partons (gluons or quark) and increasing 
its virtual mass-squared, as it approaches the hard process. In order to be plain we 
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can consider an outgoing quark from the hard process, having virtual mass-squared 
—t n , carrying a momentum fraction x n , which evolves to less virtual masses and 
momentum fraction, by emitting multiple gluons at small angles. 
The total rate of the hard process will depend on the momentum fraction distribution 
of the partons fi(x, Q 2 ) inside the hadron, as seen by an external probe (for example 
the virtual photon in the deep inelastic scattering). 

GLAPD (Gribov-Lipatov-Altarelli-Parisi-Dokshitzer) |2*5^ l2*o] evolution equations 
describe the evolution of the distribution fi(x,Q 2 ) changing the scale of the hard 
process: they are a set of coupled integro-differential equations 



t^Mx, t) = f ^PMMx/z, t). 



(2.65) 



The index j runs over the different types of partons, that is quark (q), antiquark (q) 
and gluon (g). 

The kernels in the equations, Pij(z) are called the parton splitting functions, taking 
into account different types of branching: pair production from a gluon P qg , gluon 
emission from a quark P qq and gluon splitting into gluons P gg . 
The charge conjugation invariance implies 
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They can expanded in perturbation theory as 



% if w + 



as 
2tt 



At the lowest order they take the look 
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The next-to-leading order kernels, P^p have been calculated in [2*7j . 
The symbol "+" denotes plus distributions, defined through the relation 



+ 



11C, - 2n 
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(2.66) 
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(2.70) 

/ 5(l-z).(2.71) 
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Let us notice that the plus distribution removes the infrared singularities for z — » 1: 
in fact the "+" distribution reproduces the summation between real and virtual 
diagrams, which assures the cancellation of singularities. Let us finally remark that 
the singularity for z — > is outside the domain of integration. 
The parton splitting function can be naively interpreted as the probability to emit 
a parton with momentum fraction z by the parent parton: the interpretation as 
probability implies that 



o 



dxP$(x) = 0, 
dxx[P^(x) + P^(x)} = 0, 



f dx x [2n f P$(x) + P$(x)] = 0. (2.73) 
Jo 

Actually the interpretation as probability is formally incorrect, because it does not 
handle carefully the infrared singularity for z — > 1. 

A more correct interpretation can be introduced, by defining the Sudakov form factor 



\(Q,Qo) = exp 



^ [ Q dQ> [ a s 



Qo 



[2.74) 



which turns out to be the probability to evolve from the scale Qo to Q without any 
branching |2*5] . 

A final observation is necessary for our purposes: as long as a single logarithmic 
accuracy is required, instead of the bare coupling one should consider the running 
coupling, evaluated at the transverse momentum squared, as suggested in 



a s -> a s [z(l - z)Q\ . (2.75) 

The effect of this redefinition of the scale where to evaluate the coupling constant, 
is necessary at next-to- leading order: in fact, expanding the running coupling one 
has 

a a {z(l - z)Q) = a s {Q) + ft log [z(l - z)\ + 0(a 3 s ). (2.76) 

At two loops, the logarithmic term, combined with the soft singularity 1/ (1—z), gives 
next-to-leading terms: since we have performed calculation with next-to-leading 
accuracy, these terms due to the running of the coupling cannot be neglected. 



2.10 Semi- inclusive Observables and Infrared Logarithms 

Another approach to study the properties of strong interactions is to define quanti- 
ties which are not fully inclusive, but select particular region of the phase space. 
Among these distributions are the so called shape variables , which permits to char- 
acterized the shape of the final hadronic jet: an example is the thrust , defined as 
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m 

T = max n (2.77) 
where n is an arbitrary an unitary vector. 

In practice one wants to define a variable X, which in our calculation will be the 
transverse momentum, and calculate the differential distribution 

d(X) - £ §, (2-78) 

or the partially (cumulative) integrated distribution 

I C x dT 

D(X) = ^- / dX'-=-, (2.79) 



To Jx dX> 

where X = X Q is the elastic point, that is the Born (lowest order) value of the 
observable. Let us consider the case where 

J- § = (2.80) 

and 

I f x dT 



, dX' = f . (2.81) 

To Jx dX> K } 

Radiative corrections spread the spectrum in general and introduced theoretical 
problems to face: in order to be calculable in perturbation theory, the considered 
quantity has to be infrared safe. Higher order corrections show infrared singularities 
for real and virtual emission separately: in infrared safe quantities these singularities 
cancel in the sum of real and virtual diagrams as stated in the sections above. 
Another way to state this property is requiring that the observable is insensitive to 
the emission of soft or collinear gluons, that is it is invariant under the branching 

Pi^Pj+fk (2.82) 

where the final momenta are collinear or one of them is very small. 
At the order 0(as) the calculation can be performed inserting in the phase space a 
kinematical constraint, describing the distribution we want to calculate, in terms of 
the variable of the system: 

r^^rj a ^ six '° mdn - (283) 

where Q represents generically the set of variables which parametrize the phase 
space. o(Q) is a function of the kinematical variables of the system and, in order 
to have an infrared safe distribution o(Q) has to vanish for soft emissions and for 
collinear emissions: 

for E gluon -> 





emission 
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In this kind of observables the cancellation of infrared singularities occurs, but large 
logarithms arise in every order of the perturbation theory. In fact the structure of 
the calculation at a finite order (fixed order), is 

2n 

D(X) = Gmn\og m \X - X \. (2.84) 

m=0 

These logarithms will be referred as infrared or large logarithms, because they be- 
come large when the variable approach the infrared region of the phase space. This 
occurs for X — > X , that is near the elastic point of the process and corresponds to 
the emission of soft or collinear gluons. 



2.10.1 General Structure for Single Gluon Emission 

As a remarkable example let us consider the rate for one real gluon emission. The 
contribution to the distribution in a generic variable u is 

1 rl 

2/ 



f R (u)= / duo dt \M\ 2 (u,t) 5(u-o(u,t)), (2.85) 
Jo Jo 

where |.M| 2 (u;,t) is the squared amplitude for the emission of one real gluon, inte- 
grated over the azimuthal angle and averaged/summed over the helicities and the 
colors of the initial/final partons. We defined the unitary variables u and t 

u =-=-*- (2.86) 

9 

1 — cos 9 

t = . (2.87) 

The function o(u,t) is the kinematical constraint which specifies the states of the 
infrared gluon which are included in the distribution and it is rescaled in such a way 
that 

0<o(u,t)<l. (2.88) 

For a collinear-safe observable 

o(u,t = 0) = (2.89) 

while for a soft safe observable we have 

o(w = 0,t) = 0. (2.90) 

An observable is therefore infrared safe when both conditions above are satisfied, 
or, in other words when it is unaffected by the emission of a soft and/or collinear 
gluon. If conditions ()2.89|) or ()2.90|) holds we have the cancellation of the two delta 
functions on the r.h.s. of equation (j2.100|) for a collinear or soft gluon respectively . 
As a consequence, the corresponding singularity is screened. The distribution can be 
therefore computed in perturbation theory in the limit of vanishing parton masses 
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and virtualities. 

We know that the amplitude shows a pole in u = (soft singularity), due to the 
vanishing mass of the gluon, and a pole in t = (collinear singularity), due to the 
vanishing mass of the quark. 
Therefore we can write [2] 

\M\\", t) = = + asm + + a8Fl{Uf t)} (2 . 91) 

out Ujt UJ t 

where the functions A\, Si(t), Ci(u), Fi(u,t) are defined as follows: 

At = M(0,0), (2.92) 

8l ( t) s m*LzMm, (2. 9 3) 

s M( U ,t)-mo,t)-M( U ,0) + mo,0)_ (2 . 95) 

From their definition it is clear that they are finite in the soft limit 

cu -> 0, (2.96) 

as well as in the collinear one 

t -> 0. (2.97) 

By integrating the amplitude in u; and t at the border of the phase space, large 
logarithms arise from the terms which are singular in (j2.91|) : 

D{X) = 1 + ^ [A x log 2 X + B x logX + «! + r(X )) . (2.98) 

The double logarithm represents the leading term, while the single logarithm is 
referred as the next-to- leading term, K\ is a constant and r(X) is a regular function 
in the whole phase space, vanishing for X — > X . A\ is strictly related to A%, 
typically A\ times some coefficient coming from the variable we are considering. B% 
is a combination of C\ and S\. 

The leading term is the well known Sudakov double logarithm [32] , whose coefficient 
is Ai and have both soft and collinear nature. The next-to-leading term is given by 
single logarithms of soft or collinear nature. F\{u,t) contains the terms which are 
not singular in the infrared limit. 

Let us recall that contains infrared divergences which cancel in the sum with 
virtual diagrams: 

f(u) = f R (u) + f v (u), (2.99) 

giving rise to "+" distributions, as defined above. f(u) is the physical distribution 
taking into account real and virtual emissions. 

f(u) at the border of the phase space, omitting the finite term Fx(u,t), which does 
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not product any large logarithm, may be written to 0(as) as: 

'■' V A, , S 1 (t) , C x {uj) 



f(u) = 5(1 — u) + as duo I dt 

Jo Jo 



uot ^ uo ^ I 



{8[u- o(oo,t)\ -S[u]}. 

(2.100) 



According to the conditions stated the integral of the distribution over the whole 
kinematical domain is one, i.e. the distribution is normalized: 

i 

f(u)du = l. (2.101) 

This means that the total rate is unchanged with respect to its value by the emission 
of an infrared gluon. 

Obviously the functions A\, Si (t) e C\ [uS) may be obtained with an explicit eval- 
uation of the Feynman diagrams for the emission of a real gluon. It is however 
possible to give a general derivation of these functions by using the properties of 
the amplitudes in the collinear and soft limits, as we will do in the two following 
paragraphes. 

2.10.2 Terms of Collinear Origin 

The functions A\ and C\ (to) may be obtained by considering the perturbative evo- 
lution of a light quark in QCD. 
Let us define 

z = l-u, (2.102) 
z is the energy fraction carried by the light quark, so that 

o{uj,t)=uj. (2.103) 

As we can see from condition (|2.89|h this observable is not collinear safe 9 : we regulate 
the collinear divergence with a small quark mass fi ^ 0. The denominator of the 
light quark propagator is therefore modified in (p 2 = /j, 2 ): 

1 11 

(2.104) 



(p + kf-fi 2 Q 2 uo(t + fi 2 /Q 2 
where we have defined as usual 

Q = 2p° , u = ^ 



The integration over the polar angle is 



dt — = I dt 



o t + £ J£ f 



1 „1 



'The condition 12,11)31 is instead soft-safe. 
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This result is the same of the massless case (// = 0), with the angular restriction: 

,2 



A* 

known as dead cone effect. 

The fragmentation function may be written as: 



(2.105) 



D(z;Q 2 ) = S(l-z) + a 8 f f f 

Ju 2 /Q 2 1 JO 

= 5(1 - z) + 



duj 



CO 



1 - z 



+ a s C^l-z) 



[5(1 - z - u) - 5{1 - z)\ 
(2.106) 



1 ^ 2 



Let us compare with the fragmentation function to the order as, containing the 
leading Altarelli-Parisi kernel: 



D( Z] Q 2 ) 



5(1- z) + 
5(1 -z) + 



Cpo>s { 1 + z 
2tt 

C F a s 



71 



1-z 

1 

1 - z 



log — + . . . 



2tt 



1 + z) 



(2.107) 
log ^.(2.108) 



The plus-distribution in the infrared regular kernel may be made explicit as: 



[l + z) + =l + z--S{l-z). 



This decomposition allows the identification 



CAl-z) 



71 



' 2tt 



'1 + z) 



implying that 



dz d(l-z)=-- A —. 

4 71 



(2.109) 

(2.110) 
(2.111) 

(2.112) 



This is the usual collinear coefficient appearing at the as order in processes involving 
light quarks (DIS, Drell-Yan, . . . ). 



2.10.3 Soft Emissions and Eikonal Approximation 

In processes involving light quarks only, single logarithms coming from soft emis- 
sions do not appear. This means that in these processes soft gluons can be emitted 
only at small angles, at least in the approximation where only next-to-leading terms 
are taken into account. It can be proved that, in such an approximation, they are 
emitted according to an angular ordering, behaviour known as coherence |31j . In 
processes involving light quarks only, coherence is violated at subleading level. 
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Instead in processes where also heavy quarks are present coherence is already vio- 
lated at the next-to-leading level: this means that soft emissions at large angle are 
possible even at the next-to-leading order. 

In general soft emissions can be treated according to the eikonal approximation, by 
considering the emission of a gluon from the hard on-shell partons which take part 
to the Born process. In such approximation the amplitude factorizes in the Born 
amplitude multiplied by the eikonal current: 

J»(k) = -9sJ2^h T f ( 2 - 113 ) 

i=i * 

Vi are the 4-velocities of the hard partons: vf = for massless, while v\ = 1 for 
massive partons. In the massless case the 4- velocity normalization is given by fixing 
a non trivial kinematical invariant. Tj are linear operators acting in the Fock space 
for the color of the hard partons: 

Ti\a . . . Q . . . c n ) = (Tj) , \d . . . c- . . . c n ), (2.114) 

where(Tj) c , c is the generator of the parton i. This space is the direct product of the 
single parton spaces: 

\c-i . . . Ci . . . c n ) — |ci) <g> |c 2 ) <g> • • • ® \c n ). (2.115) 

The color conservation in the QCD interactions may be written in such a scheme 
as: 

n 

^TT l = 0. (2.116) 

i=i 

where the sum is over initial/final parton colors. 

In general the notation is analogous to the one used for the angular momentum in 
quantum mechanics. 

The matrices Tj for different practical cases takes the look 

quark -> {T k ) i3 = t% 
antiquark — > (T k )ij = —t^ 

gluon -> (Tjij = if bac , (2.117) 

where f a bc are structure constant of SU(3) and t a matrices of the fundamental 
representation. 

Let us note that, in the soft case, factorization involves the momenta and the colors 
of all the hard partons. An important property of the eikonal current is that it is 
conserved: 

n 

k fi r{k) = J2T i = 0. (2.118) 
i=i 

The sum over polarizations may then be substituted with the diagonal part: 



(2.119) 
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The rate is proportional to the square of the eikonal current 



l,n 



i<j 



2vi ■ Vj m m 

1 J rp ^ rp 

Vj ■ k,Vn ■ k 1 3 



9 2 s 



E 



ky 



-T 



(2.120) 



The second sum in the r.h.s. in the equation (j2.12()j) gets contributions only from 
massive partons and therefore disappears in the well- know processes as DIS, DY, 
photon fragmentation, etc., which involve only massless partons. If heavy quarks 
are involved this term does not disappear, giving rise to logarithms of soft nature. 
In the decay of a heavy flavour into a light flavour, as we will consider, the eikonal 
current may be written as: 



J"(k) 



v ■ k 



Tq- 



n' J 



T . 

n ■ k 9 



(2.121) 



We have that v 2 = 1, n 2 = 0. We normalize the vector n on the light cone imposing 
n • v = 1, in such a way that n = 1 in the rest frame of the heavy flavour. 
In the case of our interest, with just two hard colored partons (the heavy and the 
light quark) and with an hard vertex conserving the color of the quark line, the color 
conservation reads Tq = T q . Since 



T'2 rri2 rp rp r < 

Q - 1 q - 1 Q ' 1 q - °* 



(2.122) 



the contribution for a soft gluon emission is: 

d 3 l 



fs(u) 



being v ■ I 
to (|2~86j) : 



sfc t 



f27r) 3 2cj 



v ■ P 



+ 



du I dt 



o.sCf 1 



2n ■ v 
v-l n ■ I 

cxsCf 1 



{5 [u - o (u, £)] -S[u]} 



71 



out 



71 UJ 



{S[u-o(u,t)] - 5 [u]} ,(2.123) 



u and n-k = 2ut. We define the rescaled gluon momentum analogously 

k 4 ' 



F = 



By comparing with the general expression ()2.1()()j) .we obtain: 

a s C F 



7T 



Si = Srit) 



71 



(2.124) 



(2.125) 



As already discussed in the general case, S\ ^ represent a new contribution with 
respect to processes involving light quarks only (DIS, Drell-Yan, . . . ). 
Let us note that the soft function S\(t) actually does not depend on the polar angle 
t: this fact is related to the choice of the reference frame, since the heavy quark at 
rest emits soft gluons isotropically; this is a violation of coherence at the next-to- 
leading order instead of at the NNL one, like in DIS or DY. 

By using the soft factorization one can verify the value of A\ as obtained before, 
by using the collinear factorization. This last term appears indeed in both schemes 
having the soft and the collinear enhancement. 
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2.11 Resummation of Large Logarithms 

As long as aslog 2 \X — Xq\ <C 1 the perturbative approach is reliable, however, 
when this term becomes comparable to 1 problems arise: logarithmically enhanced 
terms have comparable sizes to every order of perturbation theory 

a n s log 2n \X -X \ -0(1). (2.126) 

Obviously in this situation the perturbative approach is not meaningful: one has to 
resum at least terms such as (|2. 126(1 . finding an improved perturbative expansion 
with a larger domain of applicability. 

The theory of resummation of large infrared logarithms has been developed in the 
last twenty years for QCD processes jSSl EIE EHl CUE EHj : the resummation of terms 
like in 1)2.126)1 represents the leading or double logarithmic order. A more re- 
fined approximation is preferable to reduce the dependance on the renormalization 
scale 10 and it consists in the resummation of single logarithmic terms of the form 
a s log™ \X — Xq\. This is the next-to-leading approximation we will require for our 
calculations. In this way the region of applicability of the improved perturbative 
expansion is enlarged to aslog|X — A | < 1, a much larger region than the one 
where the fixed order calculation is valid. 

The resummation of large logarithms has been widely studied and applied to QCD 
processes for those observables which exponentiate. They have particular proper- 
ties: their matrix elements can be factorized by expressing the emission of n infrared 
(soft or collinear) gluons as the product of n single gluon emissions; moreover the 
phase space and, in particular, the kinematical constraint have to factorize in the 
same way. 

Under such conditions the perturbative expansion in the infrared region give rise to 

an exponential series, which resums logarithms with the required accuracy. 

The result of the resummation of large logarithms is accomplished by the formula 

mm 

D{X) = K(asMa s ; X) + R(X; a s ). (2.127) 

The functions involved in the formula will be the object of the calculations we will 
perform in Chapters |S] and for the process we are going to consider. Here let us 
just sum up their meaning and their role: 

• £(X; as) resums large logarithms in exponentiated form. It has a perturbative 
expansion in the form 

oo n+1 

log£(X; a s ) = J2J2 G ™ a s Lm 

n=l m=l 

(2.128) 

where L = log A. The functions ^ have the form 

oo 

gi(z) = ^2g itn z n , (2.129) 

n=0 

10 A strong dependance on the renormalization scale is in general a signal that neglected higher order terms are 
large and important. 
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so that the calculation of g\ resums leading logarithms, that is terms in the form 
OgL n+1 , g 2 resums next-to- leading contributions as a^L n , g% next-to-next-to- 
leading and so on. It is worth noting that, even if the fixed order calculation 
shows at most two logarithms for each power of as, in the exponentiated for- 
mula they are rearranged in such a way that there are only at most (n + 1) 
logarithms for n power of as- as) is a universal function, that is does not 

depend on the specific process, but only on general properties of the theory. 

• K(as), the coefficient function, is a short distance process dependent function. 
It takes into account constant terms, arising to every order of the perturbative 
expansion, which do not exponentiate. It can be calculated order by order in 
perturbation theory: 

K{a s ) = 1 + ^ia s + K 2 a s + 0(a 3 s ) (2.130) 

• R{X; as), the remainder function, is a process dependent function which takes 
into account hard contributions without any logarithmic enhancement. It de- 
pends on the specific process, the distribution, the kinematics and so on. It 
can be calculated in perturbation theory: 

R(X; a s ) = n(X)a s + r 2 (X)a 2 s + 0(a 3 s ) (2.131) 

The property 

R{X; a s ) -> for X -> X (2.132) 
holds for the remainder function. 

For our calculations we will attain next-to-leading accuracy. In order to achieve this 
level of approximation several terms need to be calculated: 

• the functions g± and g 2 to resum logarithms with next-to-leading accuracy. In 
order to reach this goal, running coupling effects should be taken into account, 
as discussed in section ()2.(ij) : 

• the value of Ki in the coefficient function has to be calculated, because the 
combination 

a S Ki ■ e Lgi{asL) ~ K X a\L 2 (2.133) 

produces next-to-leading contributions. 

Moreover, the calculation of the first term of the remainder function can be relevant 
to describe hard contributions: the remainder function is negligible near the bor- 
der of the phase space when logarithms become large, but can be relevant for hard 
emission, in the opposite limit of the phase space. 

Finally let us notice that not every observable exponentiates: however for our spe- 
cific case the conditions for the exponentiation hold and the resummation can be 
performed. In particular these requirements are the factorization of matrix elements 
for soft and collinear emissions and the factorization of the phase space. The formers 
are treated in the next section, because they involve general properties of QCD, the 
latter, instead, strongly depends on the specific observable one is considering and 
therefore will be treated in Chapter HJ 
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2.11.1 Resummation of Collinear Emissions 



The exponentiation of collinear emissions is quite direct to demonstrate, because it 
derives from general properties of the theory and it passes through the solution of 
Altarelli-Parisi evolution equations, already introduced in section (|2.9j) . 
Evolution equations can be written in compact form as in ()2.65|) 



«£/«<*.«> = Ejf 



1 dz as 



P ij (z)f j (x/z,t). 



They are difficult to solve because are integro-differential equations and can be 
reduced to simple first order equations by introducing the Mellin transform 11 : 



N 



dz z 



N-l 



In this way the convolution integral becomes a simple product 

> 

3=<}i,qi,g 



d log fi 2 



(2.134) 



(2.135) 



where 7ij,jy are the moments of Altarelli-Parisi kernels, usually referred as anomalous 
dimensions 



lij,N{as{l^)) 



dx x Pij{as{jJi )). 



(2.136) 



Now, the only residual complexity is that the differential equations ()2.135|) are cou- 
pled: the simplest case to consider is the flavour non-singlet distribution, defined 

as 

f N s = J2(U-h). (2.137) 



In this case (|2.135J) becomes 

dfNS,N{^ 2 ) 



l qq ,N {oc s (/U 2 ) ) Ins, n {^ 2 ) ■ 



d log /i 2 

The lowest order approximation for the anomalous dimension 7g g ,v reads: 

N 



(2.138) 



1 

2 + ~N(N + 



I 2T 1 - 



(2.139) 



For sake of completeness let us recall the lowest order anomalous dimensions for the 
other kernels: 



'qg,N 



T, 



R 



2 + N + N 2 
N(N + l)(N + 2) 



11 For a detailed discussion sec for example [1511251 . 
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igg,N 



2C A 



1 

12 



+ 



AT(iV-l) (AT + l)(iV + 2) ^A; 



* l' 



2 + iV + iV 2 ' 
_ iV(iV 2 - 1) 

The result of equation ()2.138|) is 



fNS,N(^ 2 ) — fNS,N{fJ>o 



a 5 (/i 2 ) 



j(0) 



where d q ° q ] N 
Moreover 



'qq.N 
27T/3 ' 



a 5 (/i 2 ) 



qq.N 



exp 




!qq,N 



q 2 1 +A)asGu ) log ^TrA, 



and if one neglect the running of the coupling 

as(<? 2 ) -> as(/^ 2 ), 
in turns out that ()2.142|) is the sum of the exponential series 



OO j 

£|(^v a s(/i 2 )log 

8=0 



(2.140) 



(2.141) 



(2.142) 



(2.143) 



(2.144) 



The insertion of the correct behaviour for the running coupling introduce additional 
logarithmic corrections. 

One can go back to the distribution f^s in the space of configurations by using an 
inverse Mellin transform 



fNsix) = ^Ic 



dN x N f NS ,N- 



(2.145) 



the integration is in general complicated and can be performed numerically. 
It is worth noting that d qqj i = 0, which corresponds to the momentum conservation, 
and dq qt N < for N > 2 which states that the non-singlet distribution function 
decreases at large x. 



2.11.2 Resummation of Soft Emissions 

The main problem related to the resummation of soft gluons in QCD is related to 
the color algebra and in general to the non abelian structure of the theory. 
Let us at first consider what happens in QED for multiple photon emission. In 
general the matrix elements for the emission of n soft photons is the Born matrix 
elements times n eikonal currents, defined in this 

J»(k) = V e (2.146) 
~ Pi - k 
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where the summation runs over the i hard emitters involved in the process as external 
legs, so that 



M\(pi, h,..., k n ) = Mi 11 J"(kj) ■ e*(%, Xj). 



(2.147) 



The rate of the process turns out to be 



a 



00 -r 



n 



d 3 ki 



I 3 



^ 1 {kjjd'iJ.vJ {.kj 



(2.148) 



the factor — is introduced to take into account the bosonic nature of the photons 
and their indistinguishability. The tensor d^ v represents the sum over polarization 
of the outgoing photons: as pointed out in section ()2.10.3|) the eikonal current is 
conserved, so that the polarization tensor can be reduced to 



d 



[IV 



9fj.v 



(2.149) 



The factors in the squared bracket of ()2.148|) are integrated over the same region of 
the phase space, so that the rate for n soft photon emission can be finally written 

as 



Or. 



I r d 3 k i n 

This shows that in QED multiple photon emissions are easily resummed into an 
exponential series, because the emission of each photon is independent from other 
emissions. 

In QCD some differences arise, due to the non abelian nature of the theory and in 
particular from the color algebra: in fact when a gluon is emitted it changes the color 
of the parent parton and consequently has an influence on the following emissions. 
Color matrices, defined in (|2.1().3j) . are included in the eikonal current, which reads 



r{k) 



Pi ■ k 



(2.151) 



Its square turns out to be 
r{k)J tM {k) = ^a s 



E T rp Pi • Pj 



(Pi ■ k) 



(2.152) 



Two kind of terms are present in the previous equation: those where the factor Tf 
appear are easy to factorize because this term is proportional to the unity matrix and 
can be treated as a c-number, so that the factorization of these terms is complete. 
This argument does not apply to the first terms in the right hand side of (|2.152j) 
because the terms Tj • Tj are in general matrices and one should evaluate their matrix 
elements between physical states: 



< M (pi)\Ti ■ TjlMoipi) > 



(2.153) 
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The problem of the exponentiation of soft emissions in QCD has not therefore a 
general solution: fortunately in some cases, among them the one we are interested 
in, this problem can be trivially solved. In particular if in the process the number of 
hard partons involved is m = 2 or m = 3 the conservation of color allows to reduce 
the scalar product between color matrices to combination of square of them, which 
are c-numbers. 

In fact if m = 2 the conservation of color reads 

Tx+T 2 = T 1 = -T 2 (2.154) 

which implies 

T\ • T 2 = Tl = T| = C F (2.155) 

where the last equality is valid if the emitters are quarks as in the case of our interest. 
For m = 3 the conservation of color simply implies 

T! + T 2 = -T3, (2.156) 

which allows to reduce the scalar products as 

Tl .T 3 = \[Tl-T}-T^. (2.157) 

For m > 3 the conservation of color implies a number of relation among the color 
matrices which is not sufficient to reduce all the scalar product to c-numbers: this 
happens since the number of independent scalar product is m ^ m ~ 1 ^ while the relations 
to reduce them are just m. For m > 3 the number of relations is less than the number 
of scalar products. 

Fortunately this complication do not appear in our calculation, because the process 
we are considering has just two hard partons. 



Chapter 3 

The radiative decay B — > X s + 7 



3.1 Phenomenology of the beauty quark 

Both from a theoretical and experimental point of view, I?-meson decays have been 
widely studied in the last years, above all after the 5-factories (Belle and BaBar) 
and CLEO began to give their results. 

The great interest that b physics arouse is due to the peculiar role of this quark 
in the Standard Model: it is the heaviest quark to hadronize, since the top quark 
decays before it can give rise to hadrons, and its physics is a very sensitive test of 
the Standard Model. It provides for example a probe to discover new physics, as 
well as a way to measure CKM matrix elements and CP violation parameters. 
In what follows we will consider decays of mesons containing a beauty quark and 
another light quark 1 , such as B = B , B + , B s , . . ., that is processes such as 

B^X + W B (3.1) 

where X represents the hadronic final states and Wb a (real or virtual) intermediate 
electroweak boson (7, W^. 

The mass of the b quark is a fundamental quantity of the Standard Model because 
it enters many processes useful to measure the parameters of the theory, such as, 
for example, CKM matrix elements. 

In spite of this importance it is still a poorly known parameter: quarks are confined 
into hadrons and they never appear as free particles, so a direct measure of their 
masses, as done for leptons, is not possible. The determination of quark masses 
is indirect and quantitative results depend on the theoretical framework chosen to 
define them. 

Recent results [H| give for the beauty running mass in the MS scheme 

m bM s = 4-0 - 4.5 GeV (3.2) 
which corresponds to the pole mass 

m b ,POLE = 4.6 - 5.1 GeV (3.3) 

1 Quarkonium systems will not be considered. 
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A complete discussion about results and theoretical problems related to the beauty 
quark mass can be found in |40j . 

Since the mass of the b quark is much larger than the scale where hadronization 
occurs, Aqcd ~ 200 MeV, 

m b > Aqcd (3.4) 

the application of perturbative QCD techniques seems to be viable. The heavy quark 

expansion technique (see for example |6lj ) allows to control non perturbative effects 

in inclusive decays and implies that the rate for hadronic decays T(B — > X + Wb) 

can be well approximated by the partonic rate Y{b — > q + Wb)- 

In this picture the spectator quark in the meson is supposed not to be involved in 

the inclusive process. Obviously in exclusive decays, where different decay channels 

are distinguished, the spectator quark has an important role, but in what follows 

we will be devoted to inclusive processes only. 

The more relevant decay of the b quark is the transition 

b^c + W* (3.5) 

(where the symbol * denotes a virtual particle) but other decays with very small 
branching ratio (~ 10 -4 — 10~ 5 ) can occur and can dominate near the border of the 
phase space, where the dominant decay is forbidden for kinematical reasons. 



3.2 Rare b decays 

Decays with very small branching ratio are referred as rare decays and are interest- 
ing because they are very sensitive to poorly known parameters of the theory, for 
example small CKM elements or test the correctness of the dynamical structure of 
the Standard Model. 

An important rare decay of the b quark is the transition 

b^u + W* (3.6) 

where the semileptonic decay 

b^u + W* ->u + l + V[ (3.7) 

seems to be the cleanest process to measure the CKM element \V U b\. 

Beside there are the so called radiative decays, rare decays with the emission of a 

photon, such as 

b^s + j (3.8) 

b^d + j (3.9) 

For example ()3.9j) is interesting because it is sensible to CP violation. Several 
theoretical problems are involved in its study, making it a challenging issue for 
particle physics. Its branching ratio is very small jH] 



6.0 10~ 6 < BR(B -»■ X d7 ) < 2.6 10~ 5 



(3.10) 
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Figure 3.1: Vertex for b — > S7 in the Standard Model. 



In the rest of this thesis we will deal with the process ()3.8|) : the world average for 
the experimental branching ratio is [HI 14*2"] 

BR(B -f X s7 ) = (3.34 ±0.38) 10" 4 (3.11) 

In general, inclusive decay rates of the h quark can be calculated in perturbation 
theory in a model independent way. However exclusive decay rates are easier to 
measure, though they cannot be calculated from first principles. 
Since the exclusive channel B — ► Kj violates angular momentum conservation, the 
simplest process to observe is B — > K*^. The first measurement of the exclusive 
rate for the final states K* (892)7 has been performed at CLEO [4*3*] . 
To measure inclusive decay rates two independent methods are used JTTJ: the clean- 
est way is to measure and sum decay rates for different exclusive channels as K'-fmc 
where n < 4 and at most one pion is a 7r°. The previous constraint are necessary to 
reduce the background but introduce some model dependence. 
The previous method can be integrated with an independent one, which consist in 
measuring of the high energy photon, emitted in the process. The photon spectrum 
has to be separated from the background of other processes and this is performed 
by using a neural network. This method was used at CLEO to perform the first 
measure of the inclusive rate for b — > 57. 

Since in the Standard Model flavour changing neutral current do not occur, a decay 
such as (13. 8J) must be mediated by a loop even at the lowest level. 
As shown in the figure 1)3. ip . the particles in the loop are a virtual W and a virtual 
u-type quark, that is up, charm or top. The different contribution of the three quarks 
in the loop is mainly due to the ratio of CKM elements: since V u b is much smaller 
than V c b and V t b the u quark contribution in the loop can be safely neglected. The t 
quark is favored by its larger mass, but c quark contribution should not be left out. 
It turns out to amount about 30%. 

The photon can be radiated by each charged line, as shown in the figure. 

Since this decay is loop mediated it is a sensitive probe of new physics: for example 
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the W boson may be substituted by a charged Higgs, in an extension of the Standard 
Model with two scalar doublets or by a supersimmetric particle. If the mass 
of new physics particles is not larger by many order of magnitude than Standard 
Model particles an effect should be detectable. However it should be pointed out 
that efforts to find signals of new physics in past years have not given results and, 
at the moment, the branching ratio for this process agrees with the prediction of 
the theory. 

Recent calculations for the inclusive branching ratio gives the prediction ^SJ HH] 

BR(B -»■ X s7 ) = (3.60 ± 0.30) 10~ 4 (3.12) 

Another useful application of b — > 37 is the measure of CKM elements: in particular 
this process may provide a way to parametrize non perturbative effects in non per- 
turbative functions that can be used as an input in ()3.7|) . allowing a precise measure 

of \v ub \ m. 

The reason why we are interested in this process is that b — > 57 shows a non trivial 
interplay between electroweak and strong interactions: the process at quark level 
is mediated by electroweak interactions, but coloured particles are responsible of 
gluon emissions which can be treated by well known tools of QCD (resummation, 
evolution equations, . . . ). 

Taking into account QCD corrections different scale appear in the process: m t ~ mw 
is the mass of the heavy (virtual) particles in the loop, mj,, the b quark mass, is the 
hard scale of the process, Aq C d is the scale where non perturbative strong effects 
appear and they can be ordered as 

m t ~ m w > m b > A QC d (3.13) 

The interplay between these scales makes the process difficult to handle: they can be 
disentangled by building an effective theory where the heavy particles are integrated 
out and an effective low energy hamiltonian is considered, like in Fermi hamiltonian 
for (3 decays. 

In the next section such an hamiltonian will be briefly described. 
3.3 Effective Hamiltonian for b — > S7 

QCD corrections to the vertex in figure (|3.1j) heavily change the lowest order result 2 
and need to be handled carefully: due to the presence of different energy scales large 
logarithms of the form 

a n s \og p {m 2 w /m 2 b ) with p < n (3.14) 

make their appearance at higher orders of perturbation theory. These logarithms are 
significant because mw 3> m b and the term — log m w jm\ is too large to be a good 
expansion parameter. Terms such as ()3.14|) can be resummed with Renormalization 

2 They enhance the rate of a factor 5: in practice this is a peculiar case where perturbative corrections are more 
relevant than the lowest order 1481 . 
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Group Equations to take them under control. 

This can be done by building an effective hamiltonian of the form 3 |l9t 150] 

AC 8 

H eff (x) = -— | V2V« Y^Cfa) dfato). (3.15) 

This way long-distance effects, both perturbative and non-perturbative, are factor- 
ized in the matrix elements of the operators Oj, while the short-distance effects are 
contained in the coefficient functions Cj (/if,), calculable in perturbation theory. 
In order to build the effective hamiltonian, the heavy particles in the loop, the W 
boson and the top quark, have to be integrated out, so that they do not appear as 
physical degrees of freedom. The integration of the two heavy particles at the same 
time corresponds to neglect strong coupling running effects between m t and mw- 
Then the large logarithms are included in the coefficient functions Cj(n) by scaling 
down the coefficients to n — Oimb). Logarithms can be resummed in the coefficients 
Cj(nb) with the required accuracy: for example the leading approximation involves 
the resummation of terms with n = p, the next-to-leading one terms with p = n — 1 
and so on. This is performed in j^l EH EU EH1 El] i n MS scheme. 
At the end, an hamiltonian describing a 5 quark theory is obtained. 
A suitable basis for the operators Oj is given by six four-quark operators,the current- 
current operators 0%, O2 and the so called penguin operators O3 — Oq, and by the 
magnetic penguin operators (electromagnetic) O-j and (chromomagnetic) 0$ j^JE0] : 



o 1 = 

6 2 = 

63 = 


{cL,p7(j,b L>a ) (s L>a 7 / ,c L>/3 ) 

(cL, a ltMb L , a ) (sl,/37mCl,/3) 


6 4 = 


q 


6 5 = 


q 


6 6 = 


q 


6 7 = 


g 

i« o^Us^SL^^bR^F^ 

l07T z 


Os = 





(3.16) 



where rribWsi^b) is the b mass in the MS scheme, evaluated at //& and q = u,d,s,c 
or b. In the operators O-j and Og a term proportional to m s is present, but we have 
omitted it because m s <C m;, and the approximation m s — > is justified. 
For several reasons a different choice for the basis can be performed and we will 



3 Small CKM matrix elements are neglected. 
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follow the definition given in This basis reads 



O x = 


(sL7^ a CL)(cL7^ a 6 L ) 


02 = 


/ — \ f — in \ 


3 = 


q 


04 = 


(s Ll ,T a b L )(J2^T a q L ) 

q 


05 = 


q 


^6 — 


q 


07 = 


g 


08 = 





(3.17) 

It is important to underline that ()3.16|) and ()3.17|) are equivalent from a physical 
point of view and can be turned one into other with proper transformations. 
The dimension of these operators is six: higher- dimension operators have coefficients 
suppressed by inverse powers of the masses of the integrated particles it and W) 
and do not contribute in first approximation. 



3.4 Lowest Order Amplitude for b — > S7 

Now to calculate the inclusive rate for b — > 57 one has to evaluate the matrix 
elements of the operators of the basis (|3.17|) between external states. 
Let us consider at first the lowest level 4 . The matrix element for a generic operator 
Oi is 

<S7|0,|6>. (3.18) 

A straightforward calculation shows that the operators 0i-06 have vanishing matrix 
elements. Moreover the operator 0§ at the lowest level describes the process b — > s g, 
so that it enters only in radiative corrections. 

At the lowest level the only operator giving contribution is O7, whose matrix element 
reads 

< S1 \6 7 \b >= mb ^u(P)^ l -±J±u{p) (3.19) 
Stt z 2 

where e is the photon polarization, the photon momentum, P^ and p M are respec- 
tively the heavy and light quark momenta. 

4 In this case the expression Born amplitude to indicate the lowest order amplitude could be misleading because, 
as we saw, the process is loop mediated even in the lowest approximation. 
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The inclusive rate is 



4C f 

To = \V*V tb \ 2 — / | < s 7 |0 7 |& > I 2 d$ (3.20) 

Let us perform the calculation in the heavy quark rest frame, where 

P" = (m 6> 5), 

for a massless strange quark. 
The inclusive rate turns out to be 

a^ Glmlml^jm,) \V tb V t *\\ 2/ 
To - ^3 C 7 M , (3.21) 

As we pointed out, the parameter rrib, appearing in physical results, should be care- 
fully considered: the lowest order rate contains both m& = m^poLE, coming from 
the phase space, and, in principle, fn b jjg(fi), coming from the effective (renormal- 
ized) operator [S3]. However, as long as QCD corrections are not concerned, the 
renormalization point is not important and we can safely choose \i = m^. As we 
will see, when strong interaction are calculated, the definition of the b mass in the 
MS scheme must be coherent with the point where the coefficient functions are 
evaluated. 



3.5 QCD Radiative Corrections for b — > s"f 

QCD corrections to the process we are considering involve the calculation of Feyn- 
man diagrams for real and virtual gluon emissions. 

Once the effective hamiltonian is built, one can calculate the n-th order QCD cor- 
rection to the rate in perturbation theory, by evaluating the contributions of gluon 
emissions, that is processes such as 

b^ s-f g 1 ...g n . (3.22) 

This can be accomplished by evaluating the coefficient functions Cj(fi) at the order 
a$ and the matrix elements 

<s 1 g 1 ...g n \6 j \b> (3.23) 

for the effective operators. To obtain sensible physical results, these terms must 
be coherently evaluated in the same renormalization scheme, for example MS, as 
usually done in literature. In particular 0:5 corrections to the inclusive rate of b —>■ 57 
are calculated in [SH] for gluon brehmsstrahlung and in [23 EH] for virtual diagrams. 
Not all the operators of the basis ()3.17|) give a contribution as relevant as the others 
to the inclusive rate: it turns out that 7 is responsible about for 85% of the total 
rate and the remaining contributions are mainly due to the operators O2 and 0%. 
In practice the other operators can be neglected because they are associated to very 
small coefficient functions or give vanishing amplitudes. One has just to pay some 
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attention to the operators O5 and (%, giving non negligible terms, which can be 
taken into account with a redefinition of the coefficients functions of C7 and Cg 

[SHUSH El 



C1...6 
c 7 



C1...6 
c 8 + c 5 . 



(3.24) 



In this way one can mainly consider the operator O7 and smaller terms given by O2 
and Og. 

The evolution of the coefficient functions (//) is governed by the renormalization 
group equation 

^Cf f (ri = lf(a s )Cf (3.25) 



where rfJJ is the anomalous dimension matrix related to the effective coefficient 

_____ 

functions defined in (|3.24j) . In recent years, efforts to achieve a next-to- leading 
accuracy has been done, to reduce the renormalization scale dependence of the 
amplitude, which was found to be the largest source of uncertainty • This request 
needs the knowledge of the (effective) coefficient functions up to 0{p.s) 



C 



eff 



C 



(O)eff 



(a0 + 



_gQ_) 

4rr 



a 



(l)e// 



(/.) + O(as) 



and the anomalous dimension matrix up to 0(a 2 s ) 

2 



eff 
r v.. 

hj 



^lJ0)eff 

4tt 1ij 



4rr 



AVeff 
1 <J 



+ 0(a] 



(3.26) 



(3.27) 



This huge work has been performed by different groups in the last years and final 
results are listed in |53*| _2] and in references therein. 

Further details about the QCD corrections are discussed in the chapters __] and 
where explicit calculations are shown. Here let us finally recall the result for the 
contribution of the leading operator 7 to the order as as calculated in [SHI __3 EH] : 



r 7 



r c?C«) 



7i 



1 + as 

7T 



16 4 2 16 m b 

7T H log — 

9 9 3 & n 

16 4 2 4, m b 

7T H — log — 

9 9 3 S n 



The two expressions in ()3.28j) reduce one another because of the relation 



b,MS 



(M) 



m 



b,MS 



(m b ) 1 + 



2a s (m b ) m b 
log — 

7T fl 



(3.28) 



(3.29) 



and can be chosen depending which mass we want in the lowest order amplitude, 
m bMs^) or n\Ms(m 6 ). 
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3.6 Photon Spectrum 

The spectrum of the photon is the first distribution calculated and measured for 
b — > 57. 

The kinematics of this process is very simple: at the lowest level it is a two body 
decay, so that the photon energy is fixed by the kinematics: in the approximation 
of a massless strange quark the spectrum is discrete and monochromatic and the 
photon energy is: 

E, = ^. (3.30) 

This approximation is obviously too naive, because it does not take into account 
two fundamental effects: 

• gluon emissions radically change the spectrum, making the process a multi- 
body decay: the spectrum becomes continuous 

< E 1 < E^ IAX = ^ (3.31) 

QCD corrections have to calculated in perturbation theory. As discussed in 
ChapterElan improved perturbative expansion is needed to have reliable results 
in the whole phase space; 

• non perturbative effects play a role: for small exchanged momenta, non pertur- 
bative effects make their appearance, related to final states effects (hadroniza- 
tion) and initial state effects (Fermi motion). Obviously they cannot be calcu- 
lated in perturbation theory and need to be factorized in some fragmentation 
or structure function. They will be discussed in Chapter |BJ They spread the 
spectrum around the endpoint. 



The photon spectrum was first calculated in [59] for a massive strange quark. A 
comparison with the data present at the time is performed in [EE] , where a model to 
deal with the Fermi motion effects is introduced. The comparison, used to extract 

the ratio 



r , agrees with the prediction of the Standard Model, based on the uni- 

tarity triangle. Latest experimental results from CLEO collaboration are shown in 
|5D] . Other details about the photon spectrum will be discussed in section f|6.1|) . 
Let us only remark that this is a distribution where large logarithms appear near the 
border of the phase space, related to emission of infrared gluons. This happens in 
the region near the endpoint of the photon spectrum, for E 1 — > These large 

logarithms spoil the perturbative expansion in that region and the usual approach 
to handle them is their resummation according to the theory described in section 

(mm. 

In section (jnj) this will be shown and the limits of this approach will be described. 
To conclude let us underline that the photon spectrum has been the first distribu- 
tion studied for the process we are dealing with: the photon spectrum is sensitive 
to longitudinal degrees of freedom, along the photon direction. On the contrary the 
transverse momentum distribution of the strange quark with respect to the photon 
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Figure 3.2: Photon spectrum in b — > S7 measured at CLEO 60 . 

direction is sensitive to transverse degrees of freedom, so that it can give comple- 
mentary information about the dynamics of the process. 

In Chapter |Hl the photon spectrum and the transverse momentum distribution will 
be compared: being sensitive to different degrees of freedom they can give comple- 
mentary informations about several features of b physics, such as non perturbative 
effects, reliability of the perturbative expansion and so on. 



Part II 



Transverse Momentum 
Distribution in b — > 57 
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Chapter 4 

Preliminary Observations 



4.1 Motivations 

In the next chapters the calculations concerning the transverse momentum distri- 
bution in the rare decay b — ► sj will be discussed and described in great detail. 
This topic was introduced in |T] and the calculation was completed in |2|: here we 
want to underline the motivations which induce us to analyze this argument and its 
relevance in the outline of b physics. 

In the last years experiments dedicated to the study of b physics have been set up: 
in particular BaBar at Fermilab and Belle at KEK began to carry out a wide pro- 
gram of measurements of parameters involved in the physics of this quark. The rich 
phenomenology of the b physics opens many possibilities both experimentally and 
theoretically: the final aim remains a precise test of the Standard Model, with the 
determination of every parameter involved in the predictions of the theory and a 
compelling comparison with experimental data. This can show the level of accuracy 
of the theory. 

Since the discovery of the b quark the importance of QCD corrections in processes 
involving it was evident: the tree level process is always mediated by an electroweak 
current, but emissions of gluon by initial or final quark change radically the theoret- 
ical prediction. The most manifest case is just related to b — > 57 which is one of the 
few processes where radiative corrections have a larger size than the lowest order. 
Perturbative QCD was widely applied to the b physics in past years and in particu- 
lar also to its decays: the b quark is the lightest quark where this kind of approach 
can be considered viable. Its mass, as discussed in the chapter |31 is larger than the 
typical scale of hadronization, see eq. (|3.4|h and the strong coupling constant as can 
be considered a good parameter for the perturbative expansion, being reasonably 
smaller than one [5 : 

a s {m b ) S 0.22. (4.1) 

For lighter quark these considerations are not valid: even the quark charm is too 
light for a reasonable application of pQCD without heavy corrections. 
This has driven to consider perturbative QCD as an essential tool in b decays: 
in particular also more refined tools of the theory have been applied to this kind 
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of processes, for example the resummation of large logarithms for semi-inclusive 
observables has become a standard technique. 

In the last twenty years QCD has had amazing confirmations of its predictions in 
experiments of high energy physics, for example at LEP: however in these cases the 
scale of energy where the processes occurred was much higher with respect to the b 
decays. Typically these experiments where performed at scales larger than Z° mass, 

M z o = 91GeV, (4.2) 

where the strong coupling constant is about an half the coupling constant at the b 
mass energy: 

a s (M z o) = 0.12. (4.3) 

Moreover in b physics theoretical prediction are complicated by the presence of 
parameters in principle unknown and that cannot be extracted by the theory: in 
particular the CKM matrix elements and the quark mass which, as discussed in 
Chapter El are poorly known. These parameters are affected by uncertainties which 
can confuse QCD predictions. They are usually eliminated by considering ratio of 
widths of different processes, which cancel the dependance on the quark mass and 
in several cases also on CKM matrix elements. A rather good theoretical quantity 
is represented, for instance, by the semileptonic branching ratio: 

Bsl = (4.4) 

1 TOT 

which turns out to be marginally in agreement with present data [61J. 
These reasons drove us to calculate with high accuracy (next-to-leading level resum- 
mation + fixed order hard contributions) a distribution in a decay of the b quark. 
The second question to answer is why we decided to study this particular decay, 
b — > 37: as stated in chapter |3] this is a rare decay, has a small branching ratio, 
presents many experimental difficulties for the measurements and therefore seems 
not to be the most immediate choice. However the reasons we have studied it is that 
has interesting and suitable features: being at the lowest level a two body decay its 
kinematics is very simple, so that one can face the problem of radiative corrections 
and resummation without additional and unimportant complications. Moreover the 
results are in many cases very general and can be extended to other processes such 
as the semileptonic decay b — > u + e + V e or the phenomenologically more important 
b — > c + e + I7 e , where the mass of the final quark is a big kinematical complication. 
The tool of resummation was mainly applied, in past years, to processes involving 
light quarks only, while in this case the presence of an heavy quark in the initial 
states gives additional and new effects. 

In our case this scheme is justified by the fact that the double logarithm appearing 
to order as can become rather large (with respect to 1 coming from the tree level): 

-^log 2 4~-0.7 (4.5) 
An mi 

if we push the transverse momentum to such small values as p t ~ ^qcd — 300 
MeV. The single logarithm can also become rather large, having a large numerical 
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coefficient: 

log — t ~ 0.6. (4.6) 

An ml 

The purpose of resumming classes of such terms therefore seems quite justified. If 
we consider running coupling effects, i.e. if the (frozen) coupling evaluated at the 
hard scale Q = m B = 5.2 GeV is replaced by the coupling evaluated at the gluon 
transverse momentum, 

a s (m b ) -> a s (p t ) = 0.45 for p t = 1 GeV, (4.7) 

the logarithmic terms have sizes of order: 

-^#^log 2 4~-0.5 (4-8) 
An ml 

and 

An ml 

The main difference with respect to resummation in Z° decays is a hard scale smaller 
by over an order of magnitude, i.e. a coupling larger by a factor 2 and infrared log- 
arithms smaller by a factor 3. 

Finally this decay, even if rare, has been widely studied experimentally because it 
seemed to be a privileged process for the discovery of new physics: matter of fact all 
the discrepancies between experimental data and theoretical predictions could be 
explained and up to now no signals of new physics have been found. However the 
large effort to measure this process have produced a huge amount of experimental 
data, which will be possibly compared with our predictions. 

The last question to make clear is why we decide to study the transverse momentum 
distribution: the main reason is that for this process the photon spectrum has been 
calculated, measured and the role of non perturbative contributions has been dis- 
cussed, with the introduction of a structure function, called shape function jHSj. The 
transverse momentum is sensitive to a complementary kinematics with respect to 
the photon spectrum, so that it can give additional informations about the process. 
In the following chapters the calculation of the transverse momentum distribution 
will be discussed: according to the theory of resummation we decided to achieve 
a next-to-leading accuracy in the calculation, for the reasons discussed in section 

dun). 

This level of accuracy requires the calculation of the coefficient controlling double 
and single logarithmic effects and the resummation of logarithms in the function 



gi and <?2 introduced in section (|2.11|) . This part, the next-to-leading-logarithmic 
approximation (NLL), is described in Chapter 

In order to achieve a complete next-to-leading order accuracy (NLO), hard terms 
to the order as have been calculated in Chapter which shows the computation of 
real and virtual Feynman diagrams to the order as- 

The resummation of the transverse momentum show different results with respect 
to the photon spectrum and a comparison between these two complementary distri- 
butions is performed in chapter El 
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4.2 Kinematics 

The specific process we are considering is 

B^h s + X + ll (4.10) 

where h s is a strange hadron, typically a K*, and X every other hadronic state. 
Without any generality loss we can consider the B-meson rest frame and define as 
reference direction the direction of flight of the photon (axis +z). 
Let us define the transverse momentum p t of the strange hadron with respect of this 
direction: it is in practice the projection of the momentum of h s on the plane x — y. 
At the lowest order of the partonic level the process is 




Pt 



Figure 4.1: Kinematics of the transverse momentum distribution. 



6^s + 7, (4.11) 

so that, being a two body decay, the photon and the strange quark are emitted along 
the same direction and the transverse momentum, as defined above, vanishes. 
However it should be reminded that the physical process that can be detected and 
measured is the hadronic one ()4.1()j) and radiative corrections change radically the 
spectrum at the partonic level, so that different sources of transverse momentum 
can be identified: 

• Multiple gluon emissions: the emission of gluons by the initial or the final quark 
makes the process a multi-body decay 

b^s igi ...g n (4.12) 

and spreads the transverse momentum spectrum. This effect can be calculated 
order by order in perturbation theory and produce a transverse momentum p p 
for the strange quark because of momentum conservation. 

• Initial state effects: a first non-perturbative effect is the well-known Fermi 
motion Classically, this is a small vibratory motion of the b quark in 
the B meson, caused by the momentum transfer with the valence quark. In 
the quantum theory, the interaction with the light quark produces also some 
virtuality of the heavy quark. Let us work in the B rest frame. We can 
parametrize the heavy quark momentum as: 



p b = m B v + p f . 



(4.13) 
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Figure 4.2: Multiple gluon emission from the initial or the final quark in the heavy quark-meson 
rest frame. 



v M is the velocity of the heavy quark, 

^ = (1; 0,0,0) (4.14) 

and pf is the residual momentum associated to the Fermi motion, such that 

\Pf \ ~ A. (4.15) 

where A is the QCD scale. We expect that Fermi motion produces a distribution 
of intrinsic transverse momenta of the b quark of order 

\p f \ ~ A. (4.16) 

Roughly speaking, we may say that transverse momenta distributions in the 
region 

M~A (4.17) 
are sensitive to the oscillations of the heavy quark in the transverse plane (x , 

y)- 

• Final state effects: a second non-perturbative effect is the fragmentation of 
the strange quark into the hadron h s . According to the idea of local parton- 
hadron duality, hadronization is expected to modify the partonic distribution 
by a quantity of the order of the QCD scale: 

W~A. (4.18) 



The total transverse momentum can be written as the sum of the perturbative and 
non-perturbative components: 



Pt = Pp + Pf+Ph- 



(4.19) 
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and the distribution is controlled by perturbative effects only as long as \p t \ Aqcd- 
Let us note that initial and final state effects involve momenta of the same order 
and in principle they cannot be separated by any energy scale. 

4.2.1 Kinematics for Single Gluon Emission 

Let us consider the kinematics at the partonic level for one gluon emission, described 
in figure (|4.5|) . 

The mass of the strange quark can be neglected. 
Let us define the fraction of energy of the gluon as 

v = —pr = — -j 4.20 
Q m b 

being the quark mass m& the hard scale of the process. 
Let us define as angular variable 

1 — cos 9 

t = , (4.21) 

where 9 is defined as the angle between the opposite of the photon direction (axis 
-z) and the strange quark direction, as shown in the figure. 

uj and t are unitary variables and the Born kinematics (vanishing transverse momen- 
tum) is achieved for u — > or t — > 0. The transverse momentum p t (since now this 
symbols will denote the modulus of the transverse momentum) of the strange quark 
can be expressed in terms of these two unitary variables: since p t = —k t , being k t 
the transverse momentum of the gluon, because of momentum conservation, it can 
be obtained from geometrical considerations that 



p\ = k\ = ml (uH(l - t)) 



(4.22) 
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Figure 4.4: Sketch of the effects of the hadronization on the transverse momentum. 




Let us define the adimensional variable 



x = 4, (4-23) 



2 ' 



so that the kinematical constraint (see section (|2.10jl ) turns out to be 

5[x- o(oj, t)] =5[x- u)h{\ - t)\ . (4.24) 

This kinematical constraint is infrared safe, according to the definition given in 
section ([2. 10)1 . since it vanishes in the soft and in the collinear limit: 

o(u, t) — > for uj — > 

o(cu,t)^0 for t->0. (4.25) 
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It turns out from kinematics that the limits on x are 

< x < J. (4.26) 



Chapter 5 



Resummation of Transverse 
Momentum Distribution 



This Chapter is dedicated to the resummation with next-to-leading order accuracy 
of the large logarithms appearing in the transverse momentum distribution near the 
border of the phase space. 

The first step is the calculation of the coefficients controlling the double and single 
logarithms. 

Then they are resummed to every order by introducing the impact parameter space 
according to the theory outlined in section (j2.11j) . 

5.1 Calculation of the Coefficients A\ and B\ 

In this section the distribution of transverse momentum of the strange quark with 
respect to the photon direction is calculated with logarithmic accuracy for one gluon 
emission. The partonic process is 

b^s+^+g (5.1) 

This computation involves the coefficients A\ and Bi which control the leading and 
next-to-leading contributions in the resummed distribution. 
As defined in section (I4.2.1j) we consider the variable : 

x=& , Q = m b (5.2) 

The distribution of interest is therefore: 

, . , 1 dT , 

d{x) = T T x (5 ' 3 > 

The kinematical constraint to impose on the radiative gluon is given in ()4.24j) 

o(u,t) = u 2 t(l - t), 



69 



70 



CHAPTER 5. RESUMMATION OF TRANSVERSE MOMENTUM DISTRIBUTION 



so that p t distribution is infrared safe. 

The aim of this section is the determination of the logarithmic structure at one loop: 
in order to extract the logarithms arising at this stage the kinematical constraint 
()4.24j) can be expressed with its lowest order approximation 

o(u,t) = ooH. (5.4) 

Using the decomposition of the matrix element discussed in section (j2.1()j) (see eq. 
(12.91)) ) we can write to the order as'- 

- 1 rl r A, | S 1 (t) | d(u) 



d(x) = 5(x) + as duo \ dt 
'o Jo 



[5( x - uH) - 5(x)] . (5.5) 

UOZ 00 X 

The first 5-function describes the real emissions, while the second one describes the 
virtual emissions. The function F x introduced in ()2.91|) can be neglected because 
it does not produce any logarithmic enhancement: it will be taken into account in 
Chapter (0). 

By integrating over the gluon energy and polar angle one obtains the following 
contributions l : 

• leading double-logarithmic term: 

as A i r^f l % X -.H)-S {x )] = -^( 1 ^-) ; (5.6) 
Jo 10 Jo 1 A \ x / + 

• next-to-leading soft term: 

a s [ dt [ du^^{5{x - uH) - 5{x)\ = ( -) + O(l); (5.7) 



JO 



00 I \ X 



• next-to-leading collinear term: 

a s duo dt^^- [S(x - uoh) - 8{x)\ = a s C x {^j + O (J-^J . (5.8) 

Let us underline that the collinear term has a different rest with respect to the 
soft one: this rest is divergent for x — > even if it does not give rise logarithmic 
enhancement 2 . 

By summing the contributions, the one soft-gluon emission distribution reads: 

/W = iW -^il(!2I£) + + asBl (I) + (5.9) 

where 

(5.io, 

In usual hard processes, such as DIS and DY, S% = so that B\ — C\. In our case, 
single logarithmic effects are more pronounced because this coefficient is almost 
factor 2 larger. 

^^Let us note that the contribution in eg. 15.61 is called double logarithmic because it becomes proportional to 
log 2 x in the integrated rate fg dx ^ . For the same reason eq. l5~7l and EHJ 

are called single logarithms. 
2 Let us underline that every contribution of the rest to D(x), eq. 12.791 is finite. 
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5.2 Exponentiation in the Impact Parameter Space 

In this section we are going to deal with the resummation of multiple soft gluon 
emission: this is performed by factorizing QCD amplitudes for multiple soft gluon 
emission and introducing the impact parameter space which factorizes the kinemat- 
ical constraint for transverse momentum [33J E2E EH OH EH] • After that one gets the 
exponentiation of one gluon distribution. 

Before dealing with the general case, let us treat in some detail the simplest non- 
trivial case, namely double emission. 

6^s + 7 + ^i + ^2 (5.11) 
According to transverse momentum conservation we have 

Pt = -hi ~ hi- (5.12) 
Multi-gluon amplitudes factorize in the infrared limit at leading level, so that 

1 d 2 T 2 11 dV^ka) 1 dT^hi) , K1 „, 

7^-77 [hi, k t 2) ^ -= -r, ^ Ti • 5 - 13 

T B dk tl dk t2 2 T B dk n T B dk t2 

where h and k 2 are gluons momenta, while Tj is the width of one gluon emission 
process. 

The distribution therefore reads: 

~ (Pt) = 5 ( Pt ) + / d Pa 5 (p t + ka) ^^7^ + (5-14) 
T dp t J T dhi 

1 / . M , L , x 1 gj (hi) 1 dT, (k t2 ) 

+- / dk tl k t2 5 (pt + k tl + k t2 ) 1 (.5.15) 

2 J T dk tl T dk t2 

The quark transverse momentum distribution is inclusive with respect to gluon 
radiation: we do not detect the transverse momenta of individual gluons, but only 
their sum, in an indirect way, by detecting the transverse momentum of the light 
quark. To solve the k t - constraint, we go to the distribution in impact parameter 
space b by taking a two-dimensional Fourier transform, 

T^^T 6 = / / dp t exp [ip t ■ b] . 5.16 

T db T dpt 

Inserting eq. (|5.14|) in the r.h.s. of eq. ()5.16|) . one obtains 



1 dT ... , 1 dT 1 1 

f- Tb {b) = 1 + r-^ (6) + 2 



1 dT1 (6) 



T db 



2 



(5.17) 



Let us now consider the general case. The transverse momentum distribution can 
be written as 

1 dT , , p , , ^ 1 dT^ , . . . 

(pt = s (pt + E f-^r- & ' 5 - 18 
r dp t ^ r dp t 
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where the distribution due to n real gluons is given by 



dp t 



n dku d V k "; h2, 'l ktn) t < / '/ 1 ' /•■ / , /.-,„). « ■ ) . .1 ! ) ) 



1=1 



dktidktz ■ ■ ■ dk 



In the infrared limit, the matrix elements for multiple emission factorize into the 
single emission ones, so that 



1 d n T n (kti, k t2 , ■ ■■k tn ) 



dkt\dkt2 ■ ■ ■ dk 



in 



1 n 
n! II 



1 dVt 
ill £j T dk u 



(k tl ) . 



(5.20) 



The factorization of the amplitudes represented by eq. ()5.20|) is a dynamical property 
of QCD. By substituting eq.flSHJ) into (l5~27IJ) and (l5~T£l) into (l5~T9T) . we obtain: 



oo q n 

n=l ' 70 1=1 



ll 5 (p t + k n +k t2 -\ h fc fn ) 

T dk t i 



1 dT(Pt 

r rfp f 

(5-21) 

Substituting eq. (|5.2()j) in eq. (|5.1()jl . we obtain the exponentiation of the effective 
one-gluon distribution in impact parameter space: 



1 dT 
T ~db 



(b) 



1 dY 1 1 
1 + ^ ( " ) + 2 



1 dV x 
T ~db 



(b) 



+ 



i.r. limit 



exp 



dT 1 
~db 



(b) 



(5.22) 



The original distribution in momentum space is recovered by an inverse transform, 

dT {p t ) 



dp t 



(5.23) 



Let us stress that the exponentiation of the one-gluon distribution is a consequence 
of two different facts: 

1. Factorization of QCD amplitudes in the infrared limit; 

2. Factorization of the kinematical constraint going to impact parameter space. 



5.3 Higher orders 

In order to achieve a next-to-leading logarithmic resummation higher order effects 
have to be taken into account. This can be performed by considering the two 
following prescriptions during the calculation: 

1. the bare coupling have to be substituted with a running coupling at two loop 
level evaluated at transverse momentum squared of the gluon [29 : 



a s -> a s {k 2 t = p t 2 ) = a s (Q 2 x), 



(5.24) 
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where the running coupling constant is evaluated with next-to-leading accuracy, 
according to equation (|2.55jl . 

= ' ' - § '°f "°f f + (PI) . (5.25) 
/3 hg{sx) \og\sx) 

s is defined as: 

2 

s ee ^. (5.26) 

The first two coefficients of the /3- function, /3 and /3 1; are independent of the 
renormalization scheme (they are universal) and have the values: 

* 33 - 2n f 3 153-19^ 

A " ~~I2^' A " 24^ ' (5 ' 27) 

where n/ = 3 is the number of active flavours. The coupling depends on the 
transverse momentum only, therefore the integration over energy and angle in 
equation ()5.5|) does not involve any running of the coupling; it is sufficient to 
make the substitution ()5.24j) in ()5.5jl . 

2. the two loop contribution in the double logarithmic term need to be included: 

A x a s -> Ata s + A 2 a 2 s . (5.28) 

where the universal coefficient A 2 is 



A 2 — — ^[Ca(— — — ) — - nt\. (5.29) 
2vr 2[ v 18 6 ; 9 /J V ; 

According to what stated in the first point one can simply substitute ()5.28|) in 

In general the resummation of transverse momentum logarithms involves the follow- 
ing two functions, both possessing power series expansion in as- 

oo 

A {a s ) = ^A n a n s = A.as + A 2 a 2 s + ■ ■ ■ , (5.30) 

n=l 

oo 

B(a s ) = Ys B n®s = Bia s + ---, (5.31) 



71=1 



5.4 Resummation of the Transverse Momentum Distribu- 
tion 

In this section we will show the results about the resummation of transverse mo- 
mentum distribution: the results shown here are based on pQ. 
In order to perform the resummation with next-to-leading accuracy we must follow 
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the prescription in section ()5.3|) : by substituting ()5.24|) and (|5.28|) equation ()5.5|) 
reads 



d(x) = 8(x) + a s !LO (Q 2 x) I duo I dt 
Si(t) , d{u) 



r A l + a s JLO (Q 2 x)A 



2 + 



+ 



L) t 



tot 

[5{x - ooH) - 5{x)] (5.32) 



After the integration over the phase space we obtain the effective distribution of 
a single gluon in the transverse momentum space, accurate to the next-to-leading 
level: 

d , s _ M \ogx | AiPi log x log log sx Bi 1 A 2 hgx 

' 2f3 x\ogsx 2j3l xlog 2 sx (3 Q x\ogsx 2(3lx\og 2 sx' 

The inclusion of virtual diagrams introduces "+" distributions. 

The resummation of transverse momentum is performed in the impact parameter 

space jSHl EH] • 

In the impact parameter space the general formula (j2.127|) becomes 

x t/r 

— — = K (a s ) S (b; a s ) + R (b; a s ) (5.34) 
1 o do 

Let us define: _ 

(5.35) 

Let us recall two important mathematical properties of ()5.35|) : it is real and it 
depends only on |fe|. 
Now we use the relation 

iOM = iij£_ i i d{x) (536) 
r dpt nT dp t 2 7rg 2 r 1 ; 1 ; 

and substitute it in the equation (j5.16|) : 

- f 1 f 2n dcj) 

d(b) = / dxd(x) / — exp liQby/x cos <p\ , (5.37) 
Jo Jo 2tt 

By combining equations (|2.ffflj) and (|5.37jl . and integrating over the azimuthal angle, 
one obtains 

d(b) = l+j dx d R (x) [J (Qb^) - 1]. (5.38) 
Jo 

Jo(s) is the Bessel function of zero order defined by the integral 

Jq(s) = / — exp [is cos0] , (5.39) 
Jo 2tt 

from which follows that Jo(0) = 1, so that the infrared singularity in x = is 
screened. 
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As shown in section (|5.2|) . by iterating for the multiple gluon emission, the single 
emission term exponentiate and we obtain 



£(&) = exp < dx dji(x) [J Q (Qby/x) — 1] 



(5.40) 



where the meaning of the function £(&), is described in section (I2.11J1 . 
In order to integrate over the x variable, the transverse momentum, we use the 
following approximation of the Bessel function jSH], such an approximation remains 
valid up to the next-to-leading level: 



where bo = 2e lE 



J (s) - 1 ~ -8(a 

1.12. We have therefore: 
-i 



£(&) ~ exp 



fc 2 /(Q 2 f> 2 ) 



da (x) dx 



(5.4i; 



(5.42) 



Where 9 (b/bo — 1/Q) in front of the integral is understood. The impact parameter 
distribution coincides with the cumulative distribution (|2.79J) in the transverse mo- 
mentum space. It is important to notice that the step approximation in equation 
(|5.41|) cuts the small transverse momenta region, since implies 



b 2 

k 2 > 
K > b 2 ' 



(5.43) 



The approximation (|5.41j) is therefore equivalent to a prescription for the non 
perturbative effects related to the infrared pole in the coupling. To order as, one 
has: 



£(&) = exp 

' 'b f 2 ,/(Q 2 fe 2 ) 

By integrating over x one obtains: 



dx -Aia s Bia s - 

2 x x 



£(&) = exp 



M , 2 
— ^-a s log 



h 2 



Bia s log 



h 2 



(5.44) 



(5.45) 



In order to obtain the resumed distribution in impact parameter to every order let 
us define dn as in eq. (|5.33|) and substitute it in the eq. ()5.42|) . Then let us integrate 
over x and write the hard scale logarithm as a function of two loop coupling, by 
using the equation 



logs 



+ 



log ((3 a s ) ■ 



(5.46) 



P a s Pi 

The impact parameter resummed distribution may be written as the exponential of 
a series of functions, as discussed in section (|2.11j) : 



£(6) = exp [Lgi(/3oa s L) + g 2 (!3 a s L) + a s g 3 (p a s L) H ] 



(5.47) 
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Figure 5.1: Plot of the function f(b) in the variable y = Q 2 b 2 /&§ (a s = 0.22) Solid line: NLO; 
dotted line: LO. 



where 



2u2 



log 



b 2 



(5.48) 



and as = «s(<5 2 )- 

We obtain the following expressions for the g\{to) and ^(w) functions: 



0iM 

02 M 



+ 



[to + log(l - to)} 
to 



(5.49) 



L 



2/?o 3 



+ log(l - w) 

+ 



1-00 
LO logfl — LO) 



+ 



^ - ~ log 2 (l - J - f log(l - c).(5.50) 

1 — LO A pi 

The expansion to order a s of the exponent reads: 

log £(&; a s ) = --A^L 2 - Bia s L - -A^o^L 3 - ^A 2 a 2 L 2 - ^/^L 2 (5.51) 
4 o 4 I 

Let us note that a single constant I?i controls the single-logarithmic effects in any 
order. The physical reason is that a soft gluon and a collinear one with the same 
transverse momenta are emitted with the same effective coupling aig(k±). 
The function £(6, as) is plotted in fig. ([5.4)1 . 

It resum the logarithmic enhancements and becomes singular when 



LO 



(5.52) 



5.4. RES UMMATION OF THE TRANSVERSE MOMENTUM DISTRIBUTION 



77 



Since 

log Q 2 b 2 

u = PoasL w bggv^' (5 ' 53) 

the singularity occurs when the transverse strange momentum becomes as small as 
the hadronic scale, 

p ± ps - ps A. (5.54) 

The singularity (|5.52|) is produced by the infrared pole (Landau pole) in the running 
coupling and signals an intrinsic limitation of resummed perturbation theory, in 
agreement with previous qualitative analysis. Let us note that the function g2 has 
basically a pole singularity in the limit ()5.52j) . while g\ has only a softer, logarithmic, 
singularity. 
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Chapter 6 



Comparison between Transverse 
Momentum Distribution and 
Threshold Distribution 

6.1 The Threshold Distribution 



For the rare decay b — > 57 the first distribution to be studied was the threshold 
distribution, that is the spectrum of the jet invariant mass and it is strictly related 
by the kinematics to the photon spectrum , so that in practice they are the same 
quantity. 

This distribution is particularly interesting both experimentally and theoretically 
near the endpoint of the spectrum, namely where 



The theoretical importance lies on the appearance of interesting effects of logarithmic 
enhancement near the endpoint, which drove to apply techniques of resummation in 
b physics. Moreover in the endpoint region non perturbative effects become relevant 
and the interplay between resummed perturbative distribution and non perturbative 
bound state effects turns out to be crucial to understand. 

From the experimental point of view physicists studying rare decays such as b — > S7 
or b — > ue~V e have to face the problem of the large background coming from the 
transition b — > c, which is overwhelming in almost the whole phase space. Since the 
mass of the quark charm cannot be neglected in the process 



m\ 



E, 



7 



2 
0. 



(6.1) 



b 



ce v e 



(6.2) 



it turns out that the endpoint of the electron spectrum is reached for 




m\ — r 
2mb 




(6.3) 
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so that in the endpoint region of the electron spectrum for b — > u 

= f (6.4) 

the dominant transition b — > c is forbidden and its background eliminated. 
An analogous problem appear for b —> sj, where the annoying background is due to 
radiative photon emissions from the tower of transitions b — > c — > s and one has to 
discriminate from strange hadrons coming from the process of interest and the ones 
coming from charm quark decays. However in the endpoint, for 

2 2 

mt — mz 

^ > -1ST' (6 ' 5) 

again the dominant transition is forbidden and the spectrum of b — > S7 can be 
cleanly measured. 

Results about fixed order calculations [SU] and resummed calculations has been 
performed for this distribution in past years and non perturbative effects has been 
included by defining a structure function, called shape function based on the 
heavy quark effective theory [51] . 

In this chapter these results will be compared with the analogous ones found for the 
transverse momentum distribution, showing peculiar differences between these two 
quantities. 



6.2 Threshold distribution at Order as 



Let us consider the energy spectrum of the photon in the rare decay b — > 57 near 
the endpoint: 



with 



m 



x 



2 • 



(6.7) 

m B m" B 

The photon spectrum then coincides with the mass distribution of the jet. The 
kinematic constraint is: 

o(uj,t)=ujt. (6.8) 
Let us note that it is infrared safe. Explicitly the distribution reads 



f(z) = 5(l-z)+a s / duo / dt 



uot uo t 



[6(1 - z - ut) - 6{1 - z)] . 

(6.9) 



By integrating over the energy and the polar angle we obtain: 



f{z) = 6(1 -z)- a s A, 



log [1 - z] 
1 - z 



+ Bxa s 



1 - z 



(6.10) 
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where the finite corrections have been neglected and we have defined 

B 1 = S 1 + C 1 . (6.11) 

Note the symmetry between the variables wet implying that the coefficient of the 
next-to-leading term is the sum of the soft and the collinear coefficient and so it 
is for io <-> t symmetries. Large infrared contributions are included by replacing 
the bare coupling with the running coupling constant evaluated at the transverse 
momentum squared 

a s -> a s k 2 , (6.12) 

where 

kl~Q 2 uj 2 t. (6.13) 
Here the symmetry between soft and collinear contributions starts to be broken. 



6.3 Resummation of the Threshold Distribution 

In this case the resummation can be performed in Mellin space: 

±rT N =f X dx x^±- f (6.14) 

r Jo r dx 

The general formula (|2.127|) for the N-th moment of the rate can be written as 

— T N = C (a s ) f N (a s ) + R N (a s ) , (6.15) 
1 o 

where now the large logarithm contains the iV-variable: 

N 

L = log — (threshold case) , (6.16) 

N 

with Nq = exp [— 7_b] . The functions g>j in the exponent are different with respect to 
the ones in the p± case and the leading and next-to-leading ones read jHSl lEH] : 

91 (A) = -^[(l-2A)log(l-2A)-2(l-A)log(l-A)], 

92 (A) = P ° A2 ~f lAl [log(l - 2A) - 2 Iog(l - A)] - ^ [log 2 (l - 2A) - 2 log 2 (l - A)] + 

+A log( i_2A) + ^l g(l-A). (6.17) 
The expansion to order a s of the exponent reads: 

log/7v = -^A ias L 2 -as {St + C x ) L-^PoalL^-^alL 2 -^ + hj^j (3 a 2 s L 2 . 

(6.18) 
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Let us comment on the above results. The single-logarithmic effects at one loop are 
controlled by the constant 

S 1 + C 1 = - 7 -^, (6.19) 

4 7T 

i.e. by the sum of the soft and the collinear coefficients, which is different from the 
p_i_-case (cf. eqs. ()5.10|) and (|5.51|) ). At two- loop they are instead controlled by a 
different constant, 

s ' + i c > = TT (620) 

The soft and the collinear terms begin to differentiate at this order and the soft one 
has a two times larger coefficient. Contrarily to the p± case, two different constants 
are needed to describe the single logarithmic effects. The dynamical difference be- 
tween soft and collinear terms is that, for a fixed jet mass, the transverse momentum 
of a soft gluon is substantially smaller that of a collinear gluon j2HllSSl 66J. 



6.4 Singularities of the Threshold Distribution 

The functions g± and g<i in (16.17}) (and therefore also the resummed distribution) 
have two different singularities [BUI IBT| f68j: 



the first one occurs when 



or, equivalently, when 



1 _ log Q 2 /m 2 x 

2 ~ logQ 2 /A 



AQ, (6.22) 

where m is the mass of the final hadronic jet s + X. In the last member of 
(j6.21|) . we have used the approximation N P3 Q 2 /m 2 x (HHI- The singularity 
1)6.21)1 signals the occurrence of non-perturbative effects in region (j6.22j) . to be 
identified with the well-known Fermi motion it is related to soft-gluon 
effects, i.e. to the terms proportional to Ai, Si and A 2 , and not to collinear 
ones, the term proportional to C\. Fermi-motion effects are therefore controlled 
by soft and not by collinear dynamics. This fact allows a factorization of Fermi- 
motion effects by means of a function taking into account soft dynamics only, 
the well-known shape function 1 [^2]. In this region, initial bound state effects 
become relevant while final-state binding effects can be neglected jEHJIH]. 



it) the second singularity occurs at 



A = 1, (6.23) 



or 



m\ ss A 2 (6.24) 

and is related to final-state hadronization effects. Both soft and collinear terms 
are singular in this region and there are non-perturbative effects related to 
initial as well as final bound-state dynamics. 



lr The shape function is also called structure function of the heavy flavours. 



6.5. THE FERMI MOTION AND THE SHAPE FUNCTION 



83 



6.5 The Fermi Motion and the Shape Function 

The main idea underlying the concept of shape function [62J is the factorization of 
initial bound state effects in a non perturbative function, based on the Heavy Quark 
Effective Theory jM], which can be extracted from a process, such as b — > 57 to be 
considered in other processes such as b — > ue~V e . 

This approach apply to the decay of an heavy hadron B, containing a b quark, into 
an inclusive hadronic state X, having a large energy and a small invariant mass; the 
kinematics of interest is therefore [7T] 



m x 



O(Aqcd) 



Ex 

Ex > Aqcd- (6.25) 

Being the heavy quark mass very large with respect to the typical scale of QCD [M] , 
the Heavy Quark Effective Theory (HQET) can be applied 

m b > E x > Aqcd- (6.26) 

The invariant mass of the final hadronic scale is small compared to its energy, but 
is large enough to justify a perturbative approach: however non perturbative effects 
related to the initial state appear and the shape function is introduced to factorize 
them. 

At the partonic level the process turns out to be 

b -> X + . . . (6.27) 

where X differ from X because it does not contain the valence quark contained in 
B. 

The heavy quark and heavy hadrons frames are the same in first approximation: 
they differ because the heavy quark has a vibratory motion inside the hadron, due 
to the exchange of soft gluons with the valence quark. 
Therefore the heavy quark momentum can be written as 

p» = m BV » + k", (6.28) 

where mg is the mass of the decaying hadron, v M is the velocity of the heavy hadron, 
taken without loss of generality in its rest frame 

^ = (1,0,0,0). (6.29) 

k 11 is the residual momentum, taking into account the exchange of gluons between 
the heavy quark and the valence quark: according to the HQET 

k^O(A QCD ). (6.30) 

The distribution of residual momentum is of non perturbative origin and repre- 
sents in practice the Fermi Motion 2 of the heavy quark inside the heavy hadron. 

2 The term Fermi Motion is taken from nuclear physics, where a similar phenomenon occurs. 
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In practice the Fermi Motion implies a change in the momentum available for the 
final state, so that the partonic invariant mass and the hadronic one differ and in 
particular 



2 

m\ = m\ + 2E x k + —^k + k 2 , (6.31) 

being k + = k$ — k z . 

Dropping terms of order 0(Aq CD ) the approximation obtained is 

m\ ~ 2E x k+. (6.32) 

This result has two main consequence: the invariant mass at partonic level is affected 
by non perturbative effects and should be calculated by a convolution with them; 
moreover the good variable describing Fermi Motion seems to be the component k + . 
Applying this formalism to the photon spectrum, for example, the distribution of 
interest is 

1 f E ,„ dT 



$(£ 7 ) = -y o <&t— (6.33) 
which can be calculated as 

$(£ 7 ) = J dk + $ pert (£ 7 - k + ) f(k + ) (6.34) 

where f(k+) is the shape function, describing the probability that the heavy quark 

has a "+" component of the momentum k + . 

According to the HQET this quantity is formally expressed as 

f(k + ) =< B\h\8{k + - iD + )h v \B >, (6.35) 

where D + = D — D z is a component of the covariant derivative, B is the B meson 

state and h v is the field describing an heavy quark with velocity v in HQET. 

An alternative way to express this quantity, avoiding 8 distributions, is to consider 

F(k + ) =< B\hl- ±—-h v \B > (6.36) 

/C_|_ — ILsjf. ~r l€ 

and the obtain the shape function by mean of the optic theorem 

f(k+) = --lmF(k + ). (6.37) 

71 

An alternative physical interpretation of the shape function is that it represent the 
probability that the decaying quark has an effective mass at disintegration time. 
In fact the momentum available for the final state is 

Q M = m B / - <f + P ~ - (6.38) 

up to terms of order 0(Aq CD ), so that one introduces the effective mass 

= m B + k + . (6.39) 
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In practice the physical interpretation is that, in order to take into account Fermi 
Motion, one can consider a decaying quark on-shell, but with an effective m* instead 
of an off-shell quark with a virtuality k + . In this way (|fi.H4j) becomes 



A final remark needs to be inserted: this approach is viable to describe effects due 
to soft gluon emissions. Hard collinear gluons cannot be described in this pattern, 
as it can be inferred by (|6.31jl . For an hard collinear gluon none of the terms in 
(j6.31|) can be neglected and the following approximation are not possible. 
In conclusion the shape function approach is based on the HQET and allows to fac- 
torize non perturbative effects due to the emission of soft gluons and it is consistent 
in the kinematical region m 2 x = O(Ex^-qcd)- It can be extracted by experimental 
data or calculated by lattice QCD [21]. From a physical point of view these effects 
are related to the Fermi Motion of the initial heavy quark. 

6.6 Comparison between Transverse Momentum and Thresh- 
old Distributions 

It is interesting the comparison between the results of section ()5.4|) for the transverse 
momentum distribution and this chapter for the threshold distribution pQ. 

• In the threshold distribution two different coefficients S± and C\ control the 
next-to-leading terms at every order, because soft and collinear logarithms dif- 
ferentiate. 

On the contrary in the transverse momentum distribution (in the next-to- 
leading logarithmic approximation) just a coefficient, indicated by Bi, is needed 
to describe single logarithmic effects. In practice soft and collinear terms act 
in the same way, with the same dynamics. 

From a physical point of view the reason that for a given invariant mass (or 
photon energy) soft gluons and hard collinear gluons have a different trans- 
verse momentum, that is a different effective coupling 3 . On the contrary in the 
transverse momentum distribution soft and collinear single logarithms cannot 
be differentiated because, tautologically, having the same transverse momen- 
tum they couple in the same way. 

• The different structure of the singularities is the most interesting comparison. 
In the threshold case two singularities are present: a first singularity, closer to 
the origin 4 is found in A = \. The soft nature is underlined by the fact that 
only the soft coefficient S\ controls this singularities (apart double logarithmic 
coefficient A 1 and A 2 ), while the collinear coefficient C\ does not appear. This 
merely soft nature suggests that it is related to initial bound state emissions, 

3 Let us recall that the effective coupling is ag(k^). 

4 Formally the region of the spectrum from E 7 = to E 1 = mi,/2 corresponds to the evolution from A = to 
A = 1 in A space. 




(6.40) 
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namely the Fermi Motion, and can be removed by mean of the shape function. 
A second singularity is found closer to the endpoint spectrum, in A = 1: this 
singularity has not a soft nature only, it is related to hadronization effects and 
cannot be faced with a shape function approach. However, in the threshold 
distribution this singularity appears after the other one, which is already fac- 
torized by the shape function and therefore is less important. 
In the p t spectrum just a singularity appears, for u — 1, which is of the same 
nature of the one discussed above for A = 1. This singularity is not related to 
Fermi Motion effects only: this reflects the considerations discussed in section 
(14. 2|) . where with naive kinematical arguments we concluded that hadroniza- 
tion and Fermi Motion effects cannot be separated in the transverse momentum 
case. 

the previous point has an immediate kinematical explanation: while in the 
threshold distribution Fermi Motion effects can be separated by a scale from 
hadronization effects (and this allows their factorization) in the p t case this 
does not occur. 

In fact the final hadronic state has a squared invariant mass [ZTj 

m 2 x =Px = (P-q) 2 , (6-41) 

being P the momentum of the heavy quark and q the momentum of the probe, 
in our case the photon. Fermi Motion effects implies a change of order A in the 
momentum of the heavy quark, so that the squared invariant mass turns out 
to be 

m x = (P ± A - qf « m\ ± 2E x k + 0(A 2 ), (6.42) 

that is a variation of order A in the momentum of the initial quark produces a 
larger effects in the invariant mass of the hadronic state 

bm\ « 2£ X A, (6.43) 

being E x = ^f. 

Hadronization effects implies a variation of order A in the momentum of the 
final hadronic state, but this leads just to a change of order A 2 in the invariant 
mass 

5m 2 x « A 2 . (6.44) 

It is evident that non perturbative effects of these two different types act at 
different scales. 

In the p t case this does not happen: a variation of order A in the initial or in 
the final state does produce a same effect, roughly speaking of the order A as 
argued in section (|4.2)1 . There is not a scale separating these two effects in this 
case. 

previous considerations imply that, since in the p t case a singularity in the 
resummed formula, having soft nature only, is not present, a shape function 
approach, introduced to factorize Fermi Motion effects, is not viable. One 
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could try a different approach based on the effective theory introduced in [72] 
or the collinear effective theory developed in |73j . but this topic has not yet 
been developed. 

• let us note that, in general, the singularities of the functions gi (u) are more 
severe than those of the functions (A) for A — > 1/2 or 1. For example, g\ [uS) 
has a logarithmic singularity, while gi (A) has an additional pre-factor 1 — 2A 
or 1 — A which softens the singularity. 

Owing to the different singularity structure, the ^-distribution is complementary 
to the threshold one and gives independent information about non-perturbative 
physics. 
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Chapter 7 

Complete 0(ag) Calculation 



This chapter is devoted to the explicit calculation of one loop Feynman diagrams 
(both real and virtual), describing the transverse momentum distribution in b — > 57. 
The results shown here are based on [2]. 

The calculation of Feynman diagrams to the order as is needed in the resummation 
program to extract the first term k\ in the perturbative expansion of the coefficient 
function K(as), introduced in equation (j2.130|) . 

This constant is necessary to achieve the next-to-leading accuracy which we have 
required, as discussed in section (|2.11j) . 

While the coefficients controlling the logarithmic enhancement can be computed by 
general properties of QCD, such as the evolution of the light quark through the 
Altarelli-Parisi kernel or the eikonal approximation, k\ is process dependent and the 
explicit evaluation of Feynman diagrams is needed. 

After the calculation of the diagrams a matching procedure must be performed to 
avoid double counting in the logarithmic terms: in practice one has to subtract the 
logarithmic terms already included and resummed in the function S(x; as), by ex- 
panding the latter to the order 0:5. 

Moreover the calculation of real diagrams can allow the calculation of the first term 
ri(x) in the expansion of the reminder function R(as] x) (see equation (j2.1Hl|) ): 
this term is not strictly required for next-to-leading resummation, but can become 
relevant in the hard region of the phase space. Let us underline that, due to com- 
putational difficulties, remainder functions are usually neglected, since they are not 
very relevant for the resummation program, and they are known for few processes. 
The calculation is performed inserting the operator O7 of the basis in equation ()3.17j) 
in the vertex and calculating corrections due to (real or virtual) gluon emissions: 
this is the leading approximation because this is the only operator having a loga- 
rithmic enhancement. Precisely the other operators shows infrared divergences that 
can be reabsorbed in the matrix element of O7, by writing an effective coefficient 
C7 according to 52 J. This implies that the contribution to the coefficient func- 
tion, for those operators, comes only from virtual diagrams, which are kinematics 
independent and therefore can be taken from literature: this is done in section (|7.3j) . 
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7.1 Complete Calculation of Real Diagrams 

Let us start with the calculation of gluon bremsstrahlung from the initial or the 
final quark. 

The three amplitudes coming from the product of these two diagrams can be depicted 
with cut diagrams, according to Cutkowski rules and the optical theorem They 
are shown in figures (j7.U7.2l7.3j) . The calculation of the diagram in figure (j7.1|) 




Figure 7.1: Interference diagram: product of diagrams describing the emission from the initial 
and the final quark. 

involves the calculation of the trace: 



Tr 



(P + m b )7^g(l + 7 5 )(p + tyyfry^l + 75) (-P - k + m b )^ f 



(P-k)(p-k) 



(7.1) 



where the operator O7 was inserted in the vertex, according to the equation (j3.19|) . 
The attribution of momenta is given in the figures. 

Let us recall that the mass of the strange quark is neglected, so that p 2 = q 2 = k 2 = 
and P 2 = m\. 

The calculation of the trace is performed in dimensional regularization, in order to 
regularize infrared divergencies, in n = 4 + e dimensions, with e > 0. 
The trace can be calculated according to the rules written in appendix |XJ): here 
n = 7 M n M . The ladder diagram in figure ()7.2|) produces the contribution 



Tr 



(P + m b )Y(P -k + m 6 )7^g(l + 7 5 )p7m^( 1 + 7b) (-P - k + m b )-f^ 



(P ■ k) 2 

Finally the diagram ()7.3j) involves the trace 



(7.2) 



Tr 



(P + m b )yq(l + 75) (p + k)Yin P {p + k)^q(l + 75) 
(p ■ k) 2 



(7.3) 
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Figure 7.3: Interference diagram. 



After the calculation of the traces 1 we can insert the kinematics of section (|4.2.1j) . 
by defining: 



t 



2£ 7 
m b 

1 — cos 6 



2 

V = ^- (7.4) 

m b 

The conservation of the quadrimomentum implies that 

»= , ty (T.5) 

1 — 0l>{± — t) 

so that the matrix element, to be integrated over the phase space, can be expressed 
in terms of just two variables, uj and t. 

Summing the real diagrams the contribution to the matrix element in dimensional 

1 Tedious details and results about this topic are not important. 
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regularization turns out to be: 

M(u,t) = (7.6) 



1 



r(l + e/2)4(l -cu(l-t)) 3 (e + 2) 



1-u 



(l-w(l-*)) 



(-8 + 8(1 - t)u - 4u 2 {l +t)- 4e(l - (1 - t)u + u 2 {l + t)) - eV(l - t)) } 
The quantity of interest is (see equation ([4.24)1 ) 

1 dr l M(U,t) j -r- I \ V O /, M 

d$(u),t) is the phase space described in appendix iBl 

The decomposition in ()2.91|) allows to underline terms giving rise to double loga- 
rithmic enhancement, (collinear or soft) logarithmic enhancement and finite terms. 

M = whj2 Y (7 - 8) 



cm 

S l (t) 



uj e (u(e + 2) -4) 
4r(l + e/2) 

l-(l-i) 1+e / 2 
T(l + e/2) ; 



„ , , Lut(2uH 3 - lu(u 3 - 2uo 2 + 8lu- 6)t 2 ) 

ri lU), t ) = ; ; — h 

V ' ; 2(1 -u{l -t)) 



+ 



ut(4u A - 10u 3 + Ibu 2 - 13w 2 + 4)i - (1 - uj) 2 (3lu 2 - 2uj + 2)) 

2(1 -u(l -t)) 



The function Fi, since it does not give rise to infrared divergencies, can be calculated 
for e — > 0. 

We are interested in the cumulative distribution, which takes from real emissions a 
contribution of the form: 



pX pi pi 1 J-p 

D R (x)= / dx' du dt — — (u,t;e) 5[x' - u 2 t{\ - t)}. (7.9) 
Jo Jo Jo r dx' 

After the integration over the phase space we expect four kinds of terms: 

• roles in the regulator e: they parametrize the infrared singularities and cancel 
in the sum with virtual diagrams because the distribution we are dealing with 
is infrared-safe; 

• Logarithmic terms diverging for x — > 0: with the matching procedure they are 
subtracted and inserted in the function £(x; as)', 

• Constant terms: they enter the coefficient function K(as); 

• Remainder functions: terms that vanish in the limit x — > 0. 
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0.2 0.4 0.6 0.8 1 



Figure 7.4: Region of the phase space to integrate: vertical axes — > oj, horizontal axes — > t. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 



Figure 7.5: Whole phase space. Figure 7.6: Complementary region. 

Integrating over x' we have: 

f 1 f 1 1 (ir 
D R (x) = du dt — — (u,t;e) 9{x -u 2 t(l -t)}. (7.10) 

Jo Jo r dx 

The remaining integrations are non-trivial because of the simultaneous presence of 
the kinematical constraint (see figure (|7.4j) and by the dimensional regularization 
parameter e. By using the identity 

6[x-uH{l -t)\ = 1 -9[ujH{1 -t) -x], (7.11) 

we separate these two effects and rewrite the distribution Dr( X^j clS db difference 
between an integral over the whole phase space (fig. 17. 5|) and a integral over the 
complementary region (fig. I7.6J) : 

f 1 f 1 1 dT f 1 f 1 1 dT 

D R (x)= / duo I dt— —(u,t;e) - / du / dt— — (to, t; 0) 9{uj 2 t(l-t)-x] + 0(e). 
Jo Jo r dx J J r dx 

(7.12) 

In this way we have separated the two main sources of difficulty in the computation 
of equation ([7. 10)1 . 
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The first integral must be evaluated for e 7^ because it contains poles in e, but is 
done over a very simple domain, independent of x: the integration can be performed 
without the introduction of special techniques with an advanced program of symbolic 
calculation. The result can be expanded in powers of e and at this stage double and 
single poles in the regulator and constant arise. 

The second integral does not contain any pole in e and therefore one can take the 
limit e — > in the integrand: it strongly depends on the kinematical constraint, 
which makes the analytical computation very complicated. 

At first let us integrate over the variable u>, since the analytical dependence on it in 
the constraint (|4.24j) seems to be simpler. The limits of integration are 

<u<l. (7.13) 



Then, after the change of variable y — t — 1, the integration can be performed over 
y, with the limits of integration 



r<y<-, (7.14) 

T 



being 



r= l ~^Z^. (7.15) 

1 + y/l-Ax 

The analytic integration is now very complicated and seems to resist to the attack 
of usual techniques: a powerful tool that can overcome this difficulty is the use of 
harmonic polylogarithms [75]. The general properties of this approach are outlined 
in appendix lEl 

The first step is the introduction of a proper basis of function: the one we singled 
out is 



1 

y 



y + l 
1 



g[0;y] 
g[-i;y] 

9[ ~ 2]V] ~ Vvi^y) 

9[ -* y] s - 2(1 (7 - 16) 

According to the theory described in appendix El the harmonic polylogarithms 
(HPL) of weight 1 are defined as: 

rv 

J[a; y] = dy' g(a; y') for a ^ 
Jo 

J[0;y] = logy. (7.17) 
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In terms of usual functions, they read: 

J[-l;y] = log(l + y) 
J[— 2;y] = 2arctan v /y 

J[-3;y] = log(l-v^v^) (7-18) 

The HPL's of weight 2 are defined for {u, v) ^ (0, 0) as 

rv ry' 
J[u,v;y]= dy' g[u;y'} dy"g[v;y"} (7.19) 
Jo Jo 

and J[0,0;y] = 1/2 log 2 y. HPL's of higher weight may be defined in an analogous 
way. They will not be used here. 

The integration is not trivial: it is performed through chains of integration by parts 
of every terms, using the property 

-^J(m w ;x) = g(a;x)J(m w _ 1 ;x) (7.20) 
and algebraic manipulations. For example 

y n J(a,w;y)dy = — —J(a,w;y)- — — g(a;y) J(w;y)dy (7.21) 
71 + 1 J n + 1 

The remaining integral can be simplified by using partial fractioning and following 
integrations by parts. 

After this process every function appearing in the result is expressed in terms of 
algebraic functions or harmonic polylogarithms: in particular every transcendental 
function is expressed combination of hpl's. 

Once the undefined integral is calculated, one has to substitute the limits of equation 
17.141 and the final result is, summing the two integrals: 

6/2 1 



D R M = C F as ' ™> 



7T \ATTj2 2 J T(l + e/2) 



2 5 1 , 2 5 . 1 
- -— -- log x- - \ogx+ - +r 1 (x) 



(7.22) 

where r\(x) is a function vanishing for x — > 0. 

From the matching procedure ri(x) can be recognized as the first term in the ex- 
pansion of the remainder funtion 2 (see section (jZ3J)) and turns out to be 

_ (r - 1)(49t 8 + 468r 7 + 1797r 6 + 3642r 5 + 4450t 4 + 3642t 3 + 1797r 2 + 468r + 49) 
ri(j) ~ 12(r + l) 5 (r 2 +3r+l) 2 

-5 - 61r - 317r 2 - 912r 3 - 1622r 4 - 1934r 5 - 1622t 6 - 912r 7 - 317r 8 - 61r 9 - 5r 10 
+ 4(t + 1) 6 (t 2 +3t + 1) 2 ° gT 

J[0, -3, r] + J[0, -3, 1/r] - 2 J[0, -1, r] + J[-l, 0, r] + J[-l, -3, r] - J[-l, -3, 1/r] 
n , „ (t + 1)(2t 2 + 7t + 2) tt _(t + 1)(2t 2 + 7t + 2) 49 

" 2 v^arctan(v^ i (T2 n +3T+1)2 (^Vsr+l) 2 +12 

+ ^logT-|log(r + l) + log 2 (r + l). (7.23) 



2 Let us remember that virtual diagrams do not depend on the kinematics, so that they cannot contribute to the 
remainder function. 
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Let us notice that r behaves as x for small values of the transverse momentum 

t(x) =x + 0(x 2 ) (7.24) 

and it is a unitary variable 

r^O for (7.25) 
r-»l for ar-^1/4. (7.26) 

(7.27) 

The relation (|7.15j) may be inverted as 

x = (7TTF- (7 ' 28) 

One can easily check that r(r) vanishes for r ~ x — > 0, by using the properties 

J[0, -1, 0] = J[0, -3, 0] = J[-l, 0, 0] = J[-l, -3, 0] = (7.29) 

7T 2 

lim J O, -3, 1/r] = lim J -1, -3, 1/r] = - — . (7.30) 

Finally let us observe that the extraction of the remainder function is based on 
a matching procedure: in order to avoid double counting we have subtracted the 
logarithmic terms already included in the resummation of logarithms, namely the 
equation ()5.9|) after the integration over x, 

— 0is '^ 1 log 2 x + atsBi log x. (7-31) 

The constant appearing in 1)7.22)1 has to be summed to the virtual contribution to 
get the coefficient function. 



7.2 Complete Calculation of Virtual Diagrams 

Virtual corrections to b — > 57 have been calculated in |qT) IoT^] for a massive strange 
quark; we present here the computation in the massless case. 

The diagrams consist of self-energy corrections to the heavy and light lines (see fig. 
(J7.7)) and 1)7.80 ) and of vertex corrections to the operator (D 7 (see figure (|7.9jl ) : we 
compute them in the MS scheme so as to be consistent with the (known) coefficient 
functions Cj. The computation can be done with standard Feynman parameter 
technique [THj or by a reduction using the integration by part identities [76J. In 
order to show both techniques, let us start by calculating the self-energies with 
Feynman parameters. The diagrams are calculated in n = 4 — e dimensions, to 
regularize ultraviolet divergencies, and off-shell to regularize infrared divergencies, 
p 2 ^ vn? . The renormalization provide for the cancellation of ultraviolet divergencies 
by counterterms in a proper scheme, decided once for all. All the calculations can 
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P 



Figure 7.7: Heavy quark self-energy diagram. Figure 7.8: Light quark self-energy diagram. 

be consistent with the chosen renormalization scheme. 
The calculation of the diagram in figure (|7.7jl gives 

d n k 1 1 



(2tt) w k 2 (P - k) 2 - m\ 



\x is the renormalization scale. 
The diagram can be written as 



Y(P-k + m b ) lfl -(Z Q -l)m b -(Z 2 -l)p. 

(7.32) 



£ = A{P 2 )+B{P 2 )P 



(7.33) 



where A(P 2 ) and B(P 2 ) can be easily calculated, by introducing Feynman param- 
eters as in appendix O They turn out to be 



A(P 2 ) 
B(P 2 ) 



47ra,c 



nm b 



dk' 

(2tt)« J \k< 2 + [P 2 - m 2 )x - P 2 x 2 ' 2 
dk' f 1 (2-n)(l-x) 



- (Z - l)m b , 
1 -(Z 2 - l)m 6 .(7.34) 



(27r)« J q [jfe/2 + (p2 _ m 2) x _ p2 x 2 

After a straightforward calculation one obtains the value of the counterterms which 
cancel the ultraviolet singularities of the diagram, in MS scheme they turn out to 
be 



(Z 
(Z 2 



I) 



as_ 

7T 

as_ 

7i 



— 7 + log 47T 



1 7 1, , 
2^ + 4-4 l0g47r 



(7.35) 



Once the ultraviolet divergencies have been cancelled, the diagrams can be calculated 
on-shell and the infrared singularities can be regularized by taking n = 4 + e. 
The contribution given by the diagram to the rate is 



Si 



OA 

dP 2 



dB 



which reads 



v as 
U\ — — 

7T 



P 2 =m 22m b + B{m 2 ) + 2m 2 

1 7 1 , ml 1 . m b 

- + -L + - log —A - 1 + t log -4 



21 



Att/j, 2 



(7.36) 
(7.37) 
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For the massless self-energies in figure (J7.8J) the case is more involved. The only term 
contributing in Si is -B(O). We need to calculate just B{p 2 ) off-shell, subtract the 
ultraviolet divergencies with the counterterm and to determine -B(O), which reads 



B(0) 



v "5 

Zji — 

7T 



1 7 1 

-2^4 + 4 l0g47r 



(7.38) 



Let us now briefly describe the evaluation of the vertex correction within the second 
method. One has to compute the scalar integral: 

d n k N(k 2 ,P ■ k,p- fc;e) 



V 



(2tt)» k 2 [(k - P) 2 - m 2 ](k - p)' 



(7.39) 



where 



N(k 2 ,P-k,p-k;e) = 32P- k-32p- k-16m 2 b + 0(e 2 )k 2 + 0(e 2 )P-kp-k + 0(e 2 )(p ■ kf . 

(7.40) 

V has at most a double pole in e coming from the product of the soft and the 
collinear singularities. The terms in the numerator N, which vanish in the soft limit 
k^ — > 0, do not give rise to soft singularities and therefore produce at most a simple 
pole coming from the collinear or the ultraviolet region. Therefore the 0(e 2 ) terms 
in iV do not contribute in the limit e — > 0. Moreover, since the only terms in N 
which can give ultraviolet divergencies are quadratic in k and can produce at most 
a single pole in e, V can be safely calculated in n = 4 + e dimensions to regularize 
just infrared divergencies. 




Figure 7.9: Vertex correction diagram. 



By expressing the scalar products in the numerator as linear combinations of the 
denominators as 



k-p = -{k 2 -{k-p) 2 ), 

k-P = l -(k 2 -{k-P) 2 + m 2 b ), 

we can reduce V to a superposition of scalar integrals of the form: 

d n k 1 



T[a, b, c] = J 



(2tt)" [k 2 ] a [(k - P) 2 - m 2 ] b [{k - p) 2 Y 



(7.41) 



(7.42) 
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with a, 5, c < 1. The above amplitudes can be related to each other by identities of 
the form [76] : 

d 

fjnh __ _ n C7 AQ^i 

difef [fc 2 ]"[(A;-P) 2 -m 2 ] f '[(fc-p) 2 ] c 1 J 

with f M = /c M ,p M ,P M . By explicitly evaluating the derivatives and re-expressing 
the scalar products using eqs. (j7.41j) . one obtains relations among amplitudes with 
shifted indices. 



P 



Figure 7.10: Massive tadpole diagram. 



By solving the above identities, one can reduce all the amplitudes to the tadpole 
(figure ITTUJ) . and obtains for the integral 3 : 



V 



16 



1 



8 + 8e —J T[0,1,0], 

e J mf 



where 



T[0,l,0] = tfp. 



m 



b 

6/2 



r(-e/2) 
l + e/2 



167r \Atth 2 

The contribution associated to the vertex correction reads 

,2 \ e/2 



71 



m 



4?r/i 2 



Fll -2 



2 2 
- + --2 

e" 1 e 



(7.44) 



(7.45) 



(7.46) 



Summing self-energies and vertex corrections, and subtracting the 1/e poles accord- 
ing to the MS scheme, one obtains for their contribution to the rate Dy 4 : 



D 



v 



71 



mi 



47T/i 2 



e/2 



r i 



+ - + 41og 
2e 



(7.47) 



We have kept the factor in front of the square bracket unexpanded to simplify the 
computation of the total rate. 

The virtual corrections to the remaining operators O^-j contain only (simple) ultra- 
violet poles in e, which are removed by renormalization; their contributions to Dy 
amount only to finite constants and log mb/nb- 

3 Such a strong reduction of 3-point function to a vacuum amplitude is possible because the only scale in the 
process is the heavy quark mass m^. Virtual corrections have indeed the lowest-order kinematics P 2 = m 2 , P ■ p = 
ml/2, p 2 = q 2 = 0. 



To factorize To one has to replace "jyjg (Mb) by m b~Ms( m b) using the formula m b -^fg(^b) = rn bWs( m b)(^ ' 



3 Cj 



-)• 
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7.3 Final Result: the Coefficient Function and the Remain- 
der Function 

Summing real and virtual contributions, the transverse momentum distribution for 
the decay b — > 37 reads, to 0(as): 



D{x) = 1 + C F — 

71 



-- log 2 X— - \ogX + Ki+ Ti(x) 



(7.48) 



As expected, the result contains a double logarithm and a single logarithm of x, a 

finite term K\ and a function ri(x) vanishing in the limit x — > 0. 

By expanding the resummed formula f)2. 127|) to order as one obtains: 

D(x) = ^1 + ^ F ^ S \ ^1 - ^-a s \og 2 x + Bia s log x^j + —^^- ri(x) 
= 1 — a^log x + Bias logs H Ki H /^(x) + C>(a 5 ). 

4 7T 7T 

(7.49) 

The matching procedure consists in the identification of the values of K\ and r\(x) 
and the computation allows to we check the values for Ax and Bi, evaluated in pQ 
using general properties of QCD radiation. 

The value of the first term in the expansion of the coefficient function is 

The value of the first term in the expansion of the remainder function is given in 
equation ([7.23]) . 

As explained in previous sections, the remaining operators O^-j contribute to D(x) 
only by finite terms fj and remainder functions. Since the constants fj come from 
virtual diagrams alone, we can quote their result from [S3 EE] and present an im- 
proved formula for the coefficient function, in analogy with [77] : 

K(as) = 1 + _g_^ rj + 7i7 log _j + __ +0(Qs) (7 . 51 ) 
where 

f» = r { 

Let us remark that only the coefficients related to the operators with i = 1, 2, 7, 8 
are relevant, because the others are multiplied by very small coefficient functions 
and can be neglected: 

1 



7.3. FINAL RESULT: THE COEFFICIENT FUNCTION AND THE REMAINDER FUNCTI0NW1 



Sr 2 



-4.092 - 12.78(0.29 - m c /m b ) 



27 



(33 - 2tt 2 



(7.53) 



The analytic expressions for the coefficient functions as well as a standard numerical 
evaluation are given in |H2]- The anomalous dimension 7^ is derived from the 



coefficient of the logarithmic term in k\. The values of 7^ are 52 : 



(o) 

in 



208 416 176 152 



243' 81 ' 81 



243 



6272 4624 32 32 
81 ' 243 'T' - 1T 



(7.54) 



The calculation of the coefficient function in (|7.51|) complete the program of next- 
to-leading resummation. 

The remainder function is in general negligible for small transverse momenta, for 
larger value of transverse momenta a small change in the fixed order prediction can 
be noticed, as shown in figure (j7.11j) . The figure shows that for larger values than 
x = 0.1 the remainder function produces a correction of the order of (9(10 — 15%) 
with respect to the only logarithmic contribution. 

Obviously this happen where logarithms are already quite small and this effect can 
be taken into account near the fully inclusive region of the spectrum, but can be 
discarded in the semi-inclusive one (small p t ). 



. 8 



0.6 



0.4 



0.2 



. 05 



) . 15 



0.2 



0.25 



Figure 7.11: Comparison between the full fixed order calculation for D(x) (solid line) and the 
logarithmic approximation (dashed line). 
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Chapter 8 

Mass Effects from the Final Quark 



8.1 Motivations 

In this chapter preliminary results are shown, concerning the structure of infrared 
logarithms in heavy flavour decays, where both the decaying quark and the final 
quark are massive. 

Example of these transitions are top — > beauty, beauty — > charm, beauty — > strange 
(with a massive strange quark). 

The main phenomenological difference with respect to transitions like beauty — > up 
or beauty — > strange (if we neglect the strange quark mass) is that in this case the 
hadronic final jet branches from a massive quark. 

In order to be plain, we deal with a specific physical distribution, the photon spec- 
trum distribution in b — > sj, but it should be clear that the results apply to other 
Q — > Q transitions with few changes in kinematics. 

In particular from a phenomenological point of view it should be very interesting to 
apply these results to the transition 

b^c. 

In fact it has a much higher rate than the rare process we are considering and 
it is widely investigated from an experimental point of view. However, being the 
kinematics more complicated, we decided to start from a simpler process, anyhow 
having a great phenomenological importance. 

In this chapter we analyze the main differences with respect to the massless case, 
due to soft gluon emissions, which are affected by the the presence of a mass in the 
final state. 

8.2 Kinematics 

Let us consider the photon spectrum in the decay 
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(8.1) 
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where X s is an hadronic jet containing a strange quark. 
At the tree level in the partonic process it is 

b -> s 7 . (8.2) 

Two practical reasons drove us to consider this specific process: first of all at the 
tree level it is a two body decay and its kinematics is very simple, moreover the 
results for this decay neglecting the strange quark mass are well known and it can 
be interesting to compare the effects due to the presence of the mass. When we will 
need such a comparison, the case with the strange mass neglected will be referred 
as the massless case. 

At the tree level the photon energy is fixed by kinematics 

E - - ~2^r- (8 - 3) 

In the following calculations m s ^ 0, but 

ml 

H = — | < 1 
m 2 b 

Radiative corrections are given by gluon emissions 

6->s70i...0„ (8.4) 

which spread and shift the photon spectrum peak. 
Let us define the invariant jet mass 

Ml = ($>) 2 

i 

where pi are parton (strange quark + gluons) momenta. 

The photon energy is related to the invariant mass by the relation 



Let us define a unitary variable X 7 as 



2m h 



X 7 = J^- 2 (8.6) 



so that Xj — > 1 is the elastic limit. 
Another useful relation is 



mi -m% 



(l-X,) = ^— 2 (8.7) 
ml — m% 

where 

m\ = M\- ml (8.8) 

If m s = 

m 2 Ml 
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and we can recognize the kinematics of the massless case. 
The correct hard scale in the process of jet evolution is 

Q = 2E X (8.9) 

as discussed in |78| . Ex is the energy of the hadronic jet and the factor 2 is intro- 
duced to simplify the calculations. 

The jet energy is related to the invariant mass by the relation 

ml + Ml , 

E x = ^ 8.10 

2m b 

Because of (18. 5j) M^-distribution and photon spectrum are strictly related: in the 
following we will calculate relations for the M^-spectrum because it involves just 
the evolution of the jet, and then we will go back again to the photon spectrum 
using (j8.5j) . 

8.3 One Soft Gluon Emission 

Let us now consider as corrections and in particular one soft real gluon emission 

b^s 19 . (8.11) 

Let us define 

J%=ri + 1>£ (8-12) 

so that p 2 x = M\. 

It is well known that the rate for one real gluon emission is affected by infrared 
singularities of soft and collinear nature, which are cancelled by virtual emissions, 
for infrared-safe variables. 

Here we are focused on the eikonal approximation, described in section (|2.1U|) . As it 
was discussed there, this approximation allows to study the structure of logarithms 
arising from soft gluon emissions and it lies on the definition of an eikonal current 

J» = gs T b -^ g s T s — — (8.13) 

Pb -Pg Ps- Pg 

where is the beauty quark momentum and Tj are color generators. 
The squared matrix element is given by 

-g^r = 4na s (2T b .T S / Pb J s - T b 2 -^- - T s 2 m ' 



(Pb ■ Pg)(Ps ■ Pg) (Pb-Pg) 2 " (Ps • Pg) 



= M bs + M bb + M ss (8.14) 

The difference with respect to the massless case is that we need to calculate all the 
three terms in the previous formula, in particular the third one is different from zero 
in this case: moreover the different kinematics produces corrections to the other 
terms. 
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The color structure is trivial because of the presence of just two coloured partons in 
the hard vertex. 

The hard vertex conserves the color, so that 



which implies 



T'2 rp2 rp rji r < 

h — Jo — i h • i s — Op. 



Because of the relation (|8.5J) we can calculate the invariant mass distribution to get 
the photon spectrum: the kinematical constraint can be derived by ()8.12j) whose 
square is 



m\ = M\- 2p x . Pg = M 2 x - 2E 9 E X I 1 



AM X 

x cos 9 



Q 2 



§.15) 



where 9 is the angle between the strange quark and the gluon. 
Let us now define 



2E n 



Q 

z = cos 9 



3l 

E x 



the constraint becomes 




[8.16) 
(8.17) 

f8.18) 



Now we have to integrate the squared matrix element over the phase space with the 
kinematical constraint 



dT 



r dM\ 



X/ SOFT 



rn 



x 



Q 2 



- -u I I 



4M 



x 



Q 2 



(8.19) 

The kinematical constraint does not introduce particular changes in the limits of 
integration of the phase space, which turn out to be 



< lu < 1, 
-1 < z < 1. 



^8.20) 



Let us consider the three terms in the eikonal current separately: the dominant 
contribution, the one giving the leading logs in the massless case, is 



M 

The integrations gives 

1 dT 

T Q dMl 



PfPs 



(Pb-Pg)(Ps-Pg) 

E X +PX 



X/bs 



Cf 1 E x , 

a s 2 lo § 

7T m A x px 



E x - Px 



[8.21] 



(8.22) 
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where px = \J E 2 X — M\ . 
Recalling that 

M\ = m^ + (l-^)(l-X 7 )) 



we have 



mi( J JL + (l-X 1 )) (8.23) 

i c?r \ i up l , . , . i 

a s — — log \fi + (1 - X 7 )] (8.24) 



T dMi; fo m2- m 2 vr 1 - X 7 



up to terms with no logarithmic enhancement. 

The argument of the logarithm signals that the most interesting feature in this case 
is the interplay between the variable (1 — X 7 ) and the parameter /x; in fact we can 
clearly discriminate three different regimes: 

• if (1— X 7 ) >/iwe fall in the massless case, getting the usual double logarithmic 
structure 

log[l-X 7 ] 
1 — X 7 ; 

• if (1 — X y ) /i, that is very near the endpoint we have a quasi-collinear 
logarithm 

1-X T 

The collinear enhancement disappears because of the mass screen from the 
massive final parton: in this case we have a single log structure; 

• if (1 — Xj) ~ [A both the terms contribute and we have an interplay between 
the massless and massive case. 

The second contribution to the current that we want to calculate is 

2 

M bb = -Anas- b -^- (8.25) 

\Pb " Pg) 



After the integration we get 



dT \ C F 1 , 

a s — — r . (8.26) 



T dMl) bh 7r m 



Recalling (jHHJ) 



This result represents the usual soft contribution as expected from the massless case: 
this is the only term not affected by the presence of a massive quark in the final 
state. 

The last term in the current is new with respect to the massless case, because it 
doesn't appear if m s = 

2 

771 

M - = - 4TOs te^F (8 ' 28) 



108 



CHAPTER 8. 



MASS EFFECTS FROM THE FINAL QUARK 




0.65 0.7 0.75 0.8 0.85 0.9 0.95 



Figure 8.1: The figure shows the different regime discussed above in the double logarithmic approx- 
imation. The red line represents the full result obtained in (|8.24|) , compared with the massless case 
(green line) and with purely massive case (blue line). The picture graphically shows the behaviours 
and the approximations discussed in the section: far from the endpoint the mass is negligible and 
the massless result is reproduced; approaching the endpoint the mass effect makes its appearance. 



The integration gives 
1 dT , 



rn 



>bb 



a s - 



x 



a s - 



ir M 2 m 2 s 
C F 1 



7r m 



m 



X, 



^8.29) 



Let us note that ()8.29|) is zero if m s = as we expected from the current. 

The result shows a pure soft term, jzjTi an d a term screened by the massive quark, 

(i-x )+ At : f° r (1 — -^7) ^ I 1 ^ ne secon d term is no longer enhanced and we can 
neglect it. 

The photon spectrum for soft emissions becomes 



1 dT 



I soft 



5(1 - X 7 ) - a s A 1 



X 



log 



;i - x 7 ) + fj, 



a s 2S 1 



X, 



+ as Si 



;i-x 7 ) + / i 



Let us note that for m 2 s — > we get the massless result, while for /i ^> (1 
soft logarithm has a coefficient twice with respect to the massless case. 
We defined 



X, 



7T 



^8.30) 
Y ) the 

(8.31) 
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Cf 

Si = - — , (8.32) 

7T 



in analogy with the massless case. 



8.4 Quasi-Collinear Emissions 

The mass of the quark in the final state provides a natural cut-off for real collinear 
emissions: gluon emissions have to respect an angular restriction, known as dead 
cone effect, as discussed in section (J2.10D . 

In this case we can tell two main behaviours: if the mass of the final quark is 
comparable to the hard scale (as maybe it could be for the transition b — > c) then the 
collinear logarithmic enhancement, matter of fact, disappears. In fact logarithmic 
terms will surely appear as in every case where different energy scales are involved, 
but these logarithms can be considered small, or to be more precise, not enhanced, 
because the two scales are not very much different. 

Otherwise, if the final quark mass does not vanish, but is small with respect to the 
hard scales, collinear logarithms turn into less divergent terms, which we will call 
quasi- collinear. 

The formal difference is that in this case the collinear logarithms does not diverge 
when the variable we are considering (for example E 7 ) reaches the endpoint of the 
spectrum; from a physical point of view this means that the collinear emission is 
screened because of the dead cone effect. 

The emission is still described by the Altarelli-Parisi kernel (I2.fi8j) . 

2 — 2uJ + UJ 2 leading 2 

Pgqi^) = -» 8-33 

UJ UJ 

but the integration over the polar angle is restricted according to the relation (|2.1()5j) : 

z>l--^. (8.34) 

The integration over the phase space gives the result in 1)8.24)1 which represents 
the leading logarithmic approximation and therefore has both collinear and soft 
enhancement. 

As stated above, while in the massless case the distribution formally diverges as 

Mi=^> (8 . 35) 

at the endpoint of the spectrum, in this case the divergence is less dramatic and 
goes as 

log// 



1-Xy 



(8.36) 



As discussed above if the mass of the initial and final quark are not much different, 
that is their masses does not different for one or more scales of magnitude, then 
this enhancement is no longer present: physically this means that if the masses are 
similar the collinear emission is so inhibited to be unrelevant. 
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8.5 Final Considerations 

When the mass of the final quark is considered a new parameter is introduced in the 
calculation and this makes the structure of the result richer even at the logarithmic 
level. Obviously we expect that a complete calculation would show an increasing 
number of regular functions, with respect to the massless case, depending on the 
kinematics and in particular from the ratio of the masses of the quarks. 
Concerning the logarithmic structure, which we are interested in, the introduction 
of the mass of the final quark introduces different regimes, which mainly depend on 
the value of this mass. 

For real emissions we have studied the behaviour of the eikonal current, showing 
that if the mass is negligible the result of the massless case are obtained: however, 
studying a specific case, for example the photon spectrum, it is shown that near 
the endpoint the presence of the mass changes the spectrum giving a less divergent 
behaviour, because the logarithm of collinear nature is screened. In fact, due to the 
presence of the mass, the collinear emission under a minimum angle is forbidden 
(dead cone effect). 

Moreover the mass introduces new terms in the eikonal current (which correctly 
disappear in the massless limit). The importance of each term strongly depend on 
the value of the mass: if the mass is small the effect is a correction to the massless 
case, if the mass is large the structure changes sizeably. In this case (see eq. (|8.3()jl ). 
for example, the double logarithmic enhancement disappears, the soft contribution 
has a coefficient twice the one of the massless case and the collinear emission can be 
neglected at this level of accuracy. 

The results discussed in the previous sections are only preliminary: they should be 
developed to study the resummation of large logarithms when they appear and they 
should apply to different phenomenologies to understand when the effects of the 
final mass can be relevant. 

In fact the presence of the mass manifest itself, as discussed, in different point of the 
spectrum: if the mass is small just at the endpoint, if the mass is large in regions 
far from the endpoint. If they manifest very near the endpoint it is very probable 
that they are overwhelmed by non perturbative effects and will be difficult to study 
and show. This implies that the application of this formalism should show different 
behaviours in the cases of the transition b —>■ c (where the final mass is large) and 
in b — > 57 (where the final mass is reasonably small). 

We expect that the calculation of virtual diagrams will have the main effect to 
introduce plus distribution, for example 



1 



1 



(8.37) 
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without changing the general structure of the result. 
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Conclusions 



The main topic of our work has been the calculation of the transverse momentum 
distribution of the strange quark with respect to the photon direction, in the rare 
decay of the beauty quark b — > sj. 

The physical motivations which lead us to study this particular process are outlined 
in chapter HJ there a discussion of the physical effects which lead to a generation of 
transverse momentum are described. In particular we focused our attention on the 
perturbative part, whose limits give informations about the structure of non per- 
turbative contributions: such a distribution is a typical case where large logarithms 
appear near the border of the phase space, for the emission of soft and collinear 
gluons. 

The appearance of logarithms always occurs in calculations where different scales 
are involved: however if the scales have very different sizes these logarithms becomes 
so large to spoil the perturbative expansion. 

In the past years a technique to resum them was developed, in order to obtain an 
improved perturbative formula, with a range of applicability in a wider region of the 
phase space: this technique was successfully applied to processes at high energies, 
namely of the order of intermediate boson mass. Our aim was an application of such 
a technique to processes at lower energies, in our case at the energy of the beauty 
quark mass. 

In this case theoretical predictions are quite involved even in inclusive quantities: 
several poorly known parameters appear, such as CKM matrix elements or the quark 
mass, so that a stringent check of QCD at this scale of energies is not trivial. Matter 
of fact there are not indications that QCD should fail at this scale of energies nor a 
very precise check of its reliability. 

The strong coupling constant seems to be small enough to allow a reliable pertur- 
bative expansions, and in distributions such as the one we studied the technique of 
resummation seems to be necessary and viable. However only a precise calculation, 
for example with next-to-leading accuracy, as we performed, compared with equally 
precise experimental data can suggest if the perturbative QCD is reliable or other 
effects play a crucial role. 

This thesis is divided in a first part, where general tools to deal with QCD correc- 
tions in b physics are outlined and a second part where the calculations we have 
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performed are explicitely described. 

After a description of the kinematical properties of the process of interest (Chapter 
HJ, in Chapter El the resummation of large logarithms, performed with the purposes 
indicated above, is shown. 

As usual even the resummed formula cannot be applied in the whole phase space, 
since near the infrared region non perturbative effects make their appearance and 
singularities appear in the resummed formula. These singularities give information 
about the non perturbative physics. 

In Chapter |U] the comparison with a complementary distribution calculated for the 
same process, the threshold distribution, is outlined: the comparison with the dif- 
ferent structure of singularities allows to conclude that in the transverse momentum 
case the non perturbative effects, above all the Fermi Motion, cannot be dealt with 
as in the threshold case. There a new non perturbative function (shape function), 
based on the HQET was introduced, here this approach is not viable, for the mo- 
tivations contained in Chapter |3 The main reason is that in the threshold case 
Fermi Motion effects and hadronization can be separated by a scale, while in the 
transverse momentum case this does not occur. One should try in the latter case 
an alternative approach, possibly based on a different effective theory, as outlined 
in the text. 

Finally the calculation is completed in Chapter where a complete computation 
of real and virtual diagrams to order as is performed: this is necessary to extract 
some ingredients appearing in the general resummed formula, above all the con- 
stants involved in the coefficient function. Moreover it allows the extraction of the 
remainder function, the regular function describing the behaviour of the distribu- 
tion in the hard region of the phase space: its knowledge is not strictly required 
for a next-to-leading resummation, but can improved the quality of the prediction 
far from the semi-inclusive region. However most of the events are expected in the 
region of small transverse momenta and therefore the phenomenological impact of 
the remainder function is expected to be of a few percent. 

On the other side, the analytical computation of this term, which could seem aca- 
demical, is an interesting application of techniques of advanced analytical calcula- 
tion, as described in Chapter [3 which can be used for future purposes to approach 
the computation of other processes. 

In the last Chapter, |HJ a different topic is sketched: the effect of the mass of the final 
quark in b decays. The result there presented are preliminary, since a resummation 
of the logarithms has not been performed yet and only the logarithmic structure 
to one loop is analyzed. However this is the first step towards a more general ap- 
proach, taking into account the logarithmic enhancement in decays with a massive 
final quark. This will have an important application for the transition b —>■ c, where 
the mass of the charm quark cannot be surely neglected. This process is very in- 
teresting from a phenomenological point of view: having a large rate, its study is 
intense at the 6-factories. 

To summarize, the main object of our work has been the application of the log- 
arithmic resummation to the decay of an heavy quark. The calculation has been 
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performed for the first time with next-to-leading accuracy, including constants and 
remainder terms, and the comparison with experimental data could be very useful 
to verify informations about the reliability of perturbative QCD in b decays, the role 
played by non perturbative physics. 

As underlined during the text, the process of interest is a rare decay, but it can be 
considered as a first step to face other processes that are very interesting from the 
phenomenological point of view. 

This work could have several and possibly interesting extensions: first of all we hope 
that a systematic phenomenological study could be carried on within an extended 
comparison with data collected at the b factories. Obviously this would be crucial 
to argue the correctness of the approaches that are used in b decays, such as pertur- 
bative QCD. 

A second point would be the extension of the shape function approach to the p t 
case: as underlined during the text, physical reasons forbid the application of this 
approach to the p t case and one should try to use different tools. Let us notice here 
that a good parametrization of non perturbative effects is crucial, not only for this 
particular decay, but also for other important decays, such as b — > u, because the 
same non perturbative functions are involved. 

A further point would be the application of these techniques to other b quark decays 
and to other distributions, for example the jet broadening, in order to reach a more 
complete knowledge of the dynamics of QCD at this scale of energies. 
Finally the results of Chapter |H1 should be completed to face more general problems, 
such as the decay of a beauty quark into another massive quark, mainly a charm, 
but also a strange. The study of mass effects in perturbative calculation and in 
the resummation of logarithms can give higher reliability of the predictions of the 
theory. 



CHAPTER 9. CONCLUSIONS 



Part III 
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Appendix A 

Feynman Rules and Numerator 
Algebra 

This appendix contains relevant Feynman rules, used in the calculations performed 
in perturbative QCD. Feynman rules are extracted by the lagrangian of the theory, 
which reads: 

Cqcd = ^(i^-m^-^Al-dvAlY + gsA^n^ 

- gf ab \d^) 2 A b »A cu -^(f eab AlA h v )(f ecd A bll A du ). (A.l) 

Ingoing quark: 



P 

Outgoing quark: 

• ► 

P 

Ingoing gluon: 

k 

Outgoing gluon: 

k 



= u(p, s) (for quark), v(p, s) (for antiquark) 

(A.2) 

= u(p, s) (for quark), v(p, s) (for antiquark) 

(A.3) 

= e(M) (A.4) 
e*(k,X) (A.5) 
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Quark propagator: 



a P (3 

Gluon propagator: 

a , f 1 k 0,v 

Quark-Gluon vertex: 



3-gluon vertex: 





a, fi 



c,p 



4-gluon vertex: 

a, \i 




-iS, 



a/3 



rPfJ, ~ m 
p2 _ m 2 _|_ ^ 



9,u + V k ~& 
k 2 + ie 



(A.6) 



(A.7) 



(A.8) 



- gsf abc [g lll/ (k-py + g up (p-qr+ 
+ 9 pfl (q-m (A.9) 



- igs[f abe f cde (g pp g U(T -g pa g up ) 

+ f ade f bce (g'"g pa -g ftf '!r] 

(A.10) 



Ghost vertices and propagator are not important for our work and will not be 
recalled. 

The factor 77 in the gluon propagator selects the gauge: all the calculation have been 
performed in the Feynman gauge which corresponds to i] = 0. 
Sum over polarizations: 



(A.ll) 



^V(p, S)v(p,s) = rfpp 



m. 



(A.12) 
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In the Feynman gauge the sum over the polarizations of external gluons allows the 
following replacement 

A 

In the calculation of matrix elements in dimensional regularization traces of gamma 
matrices arise. The following properties have been used (n = dimension of the 
space-time): 



n 

YYl, = (2-n)Y 



Tr ( 7 ^) 



Ag up + (n- A)YY- 



(A. 14) 



In n dimension one can consistently introduce a matrix 75 in a scheme where the 
property 

holds. 75 in 4 dimension reads: 



75 = «7o7i7273- 

Traces of an even number of gamma matrices and a 75 vanish: 

Tr[7s7 Ml ...7 W J = 0. 
The metric chosen to perform the calculations is 



(A.15) 
(A.16) 

(A.17) 



/ 1 











\ 





-1 
















-1 







V° 








-1 


/ 



(A.18) 



APPENDIX A. FEYNMAN RULES AND NUMERATOR ALGEBRA 



Appendix B 



Phase Space in Dimensional 
Regularization 

The formula used to calculate the decay rate of a particle in its rest frame is 

*T = ± d* n 6(£tf) |M| 2 . (B.l) 

i 

m is the mass of the particle, $ n is the n-body phase space, S(J2iPi) assures the 
momentum conservation and \M\ 2 is the squared matrix element of the process. 
For the distribution of the transverse momentum in b — > sj, \M\ 2 has been calculated 
in Chapter [7| , in this appendix the results for the 2-body and 3-body phase space 
in dimensional regularization are shown: these terms appears in the calculations of 
the Born amplitude and in radiative corrections. 
The lowest order process is b —>■ 57. 
The general expression for the phase space is 

d$2 = { f* n ( ^ n (2 7 r)^ + (p 2 )(2 7 r) ( 5 + (g 2 )(2 7 r)^(P^ - p" - <f), (B.2) 

where p M is the strange quark momentum, the photon momentum and the 
beauty quark momentum and 

S + (p 2 ) = 6(p 2 )9(p°) (B.3) 

n is the dimension of the space-time: for the regularization of infrared divergent 
digrams n = 4 + e. 

()B.2jl can be integrated over all the variable on which the matrix element do not 
depend explicitely: for example the matrix element for b — > 57 can be rearranged to 
depend only on the photon energy and all the other variables can be integrated 
out. 

At first the integration over the strange quark momentum can be performed, so that 
(IB. 2)1 turns out to be 

d^ = j^- 2 5 + [(P- q r]5 + ( g 2 ). (B.4) 
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and 



Let us recall that 

d n q K(q 2 ) = J d n q 5 + {q 2 ) = (B.5) 

d n ~ x a 1 

^ry = 2^7 _3 sin ™" 3 ^ sin "" 4 ^ ■ ■ ■ sin 9 n - 2 d9 1 d9 2 . . . d9 n . 2 dE T (B.6) 

0j are the angles with respect of the axes in the n — 1 dimensional space, for n = 4 
the formulae reduce to the well known one in physical phase space. Since the matrix 
element does not depend on the angles 9i they can be integrated, using 

Using the properties of the gamma function in appendix (0, finally 
f d n 1 q _ „_ 3 (n.-2)/2 r(n/2 - 1) _ 3 

,/ Te7" 2 * r(n - 2) ^ ^ (R8) 



Substituting (IB3j) and (IR8I) in (1B~4I) finally one gets 

rf$ 2 = J r y Zil y m »-*«J(l - y ), (B.9) 

2™- 1 W2-i r(n-2) y v y; v ; 

where the adimensional variable ?/ = 2E 1 /m has been defined and the kinematical 
property (P — q) 2 = m 2 (l — y) has been used. 

The 3-body phase space can be calculated in a similar way. It appears in the 
calculation of QCD corrections and corresponds to the process b — > S7#. 
In general it takes the look 

r]TL sJYl sJfl 1 

d ^ = (27T)" (2*)" (27r)n (27r)^ + (p 2 )(2 7 r)^(g 2 )(2 7 r)^(fe 2 )(2 7 r)^(P^ -p» - <f - 

(B.10) 

where fc M is the gluon momentum. 

The matrix element turns out to depend only on three variables: for example the 
gluon energy, the photon energy and the angle between the photon and the gluon 
and all the other variables can be integrated out. 

After the integration over the strange quark momentum p^, by using the momentum 
conservation and ()B.5j) . (jB.lOJ) reads 

d n ~ 1 a d n ~ l k 

d$ 3 = - — - — A—- — - — ——5 \ m 2 - 2m£ 7 - 2mE Q + 2E 1 E Q (\ - cos 9)] , 
(27r) n ~ 1 2E 1 (2TT) n - 1 2E g 1 7 9 7 9K n ' 

(B.ll) 

being 9 the angle between the gluon and the photon. 
Moreover 

An-\ rl n ~ 1 b on-3 n-2 

(27r)"- 1 2E 7 (2Tr) n - 1 2E g T(n-2) 7 19 9 v ' 9 7 ' 

(B.12) 
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where the angular variable t is defined as t = (1 + cos 9)/2. 

Substituting (|B.12j) in (jB.lljl and defining the adimensional energy fractions 

2E 1 2E g 

V = 1 , z = — y -, (B.13) 

m m 

the final result for the 3-body phase space turns out to be 

= 2 2n 7r n m i r(n _ 2) ^ ^ (l~t)W 6 [1 - V - Z + V z(l ~ t)} dz d V dt. 

(B.14) 
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Appendix C 



Loop Integrals in Dimensional 
Regularization 



A standard technique for the calculation of loop integrals consists in the introduction 
of Feynman parameters. This technique has been used to calculate the self energy 
corrections in virtual diagrams. 
It lies on the following identity: 

- r(n) / dx\ [ dx2 [ dx n ^ ^ x . (C.l) 



In Feynman gauge at one loop it can be used to parametrize the denominators: 



1 



dx / dy 



abc J J [ax + by + c(l - x - y)] c 



i r 1 , 1 



dx- 7T7t ^. (C.2) 



ab Jo [ax + b{l - x)] 
In Landau gauge other denominators can appear, such as: 



4r = 2 f dx 



a 2 b Jo [ax + 6(1 - x)] 

1 _ f l J x 2 



3 



a 3 6 Jo [ax + b(l - x)] 

1 12 I dx C dy ^ j. (C.3) 



a 2& c 7 Jo [axy + 6x(l - j/) + c(l - x)f 

Let us recall that Feynman gauge corresponds to choose rj = in (|A.7|) and Landau 
gauge to T] — 1. 

These parametrization are useful to calculate loop integrals in every regularization: 
in particular in dimensional regularization, after the introduction of Feynman pa- 
rameters, the integral over the loop momentum can be performed using 

d n k 1 (-1)° r(q-n/2) 1 



(2tt)» [A; 2 -C] a (16tt 2 )«/4 r(a) C°-™/ 2 
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where n is the dimension of the space-time and fc M the loop momentum to inte- 
grate. C is a constant function of the external momenta in the loop: usually the 
denominator does not appear as in (|C.4j) . but with scalar product of with ex- 
ternal momenta (k 2 + 2k ■ q + A). Such linear terms can be eliminated by adding 
and subtracting terms which make the denominator a square and by redefining the 
integration momentum with a translation: 

(k 2 + 2k-q + X) = (k 2 + 2k-q + X + q 2 -q 2 ) = ([k + qf - [q 2 - A]) -> (k'-C). (C.5) 

Another kind of integral appearing in ultraviolet divergent integrals is: 

d n k k 2 _. (-1)° nr(a-l-n/2) 1 
(XT [k 2 -C] a ~ Z (16tt 2 )™/ 4 2 T{a) ' ( ^ 

Integrals with linear terms in fc M in the numerator vanish once the reduction in (jC5|) 
is performed. 



Appendix D 

Gamma Function and 
Dilogarithms 



In the analytical evaluation of processes involving QCD corrections a wide range of 
functions appears: multi-loop calculations are very difficult to perform in analytical 
form and up to now people try to face complete calculations at two loops. 
For one loop calculations the range of function which can be encountered is well 
known: apart regular and elementary functions (powers of a variable basically), 
trascendental functions may appears, above all in semi-inclusive calculations. 
Beside logarithms, which are widely treated during the thesis and are relevant in 
the endpoint region of the spectrum, other special functions appear: in particular V 
function and dilogarithms. 

The T function basically is involved in the measure of integrals in dimensional reg- 
ularization and is defined as: 



T(x) 



y*- 1 dy, 



x > 0. 



It provides a generalization of the factorial for non integer number, since 

r(n + 1) = nT(n) = n! 
which can be extended for real positive numbers as 

T(x + 1) = xT(x). 



(D.l) 



(D.2) 



(D.3) 



In dimensional regularization T functions arise from one loop integrations and may 

give rise to poles in the regulator e. 

For example it is usual to perform the expansion 

T(e) = Z- fy + Z[^ + ^) e + 0(e 2 ), (D.4) 



1 1 



valid for e — ► 0, in case after the application of (|D.3J) . 
7 is the Euler constant , which satisfies the property 



7 



lim 

n— »oo 



11 1 

- logn+1 + - + - + ... + - 
2 3 n 



0.5772 . . . 



(D.5) 
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and 



Finally it turns out that 



7 = — / e x \ogx dx. (D.6) 



dT{x) 
dx ^ 

In one loop calculations (for example in virtual corrections) often the following 
integral appears: 

/o r (" + P) 

A useful property of the Y function is 

r(x)r(i -x) = -r^— , (D.9) 
smnx 



I 

Jo 



from which it follows 



r(i) = i 

r(i) = 0F. (D.10) 

Finally let us notice that the T function can be extended for x < 0, but in this region 
singularities appears for x G Z~ , that is for negative integer numbers (and for x — 0). 

While the T function appear in the measure of integrals in dimensional regulariza- 
tion, other special functions appears in the final analytical results: they are diloga- 
rithms. 

The dilogarithm is defined as: 

oo 2 

Li 2 (x) =J2^ for N < L ( D - U ) 

n=l " 

and has a branch cut discontinuity in the complex plane from 1 to oo. 
Equivalent definitions can be stated in integral form: 

Jo t 
Li 2 (x) = - / -rft 

1 — I 

Ll * w = i'r 2 !^*- (D12) 
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In figure (|D.1|) and ()D.2j) the real and imaginary part of the dilogarithm are depicted. 
The dilogarithm has the properties: 



Li 2 (0) = 
Li 2 (l) = 



TV 2 



6 

2 



Li 2 (-1) 



7T 

12 

' 2 , 



vr 2 1 



Li 2 (l/2) = ---log 2 2. (D.13) 



Other relevant properties are: 

7T 2 1 

Li 2 (— x) + Li 2 (— 1/x) = — - — -log 2 a; x > 

o 2 

7T 2 1 

Li 2 (x) + Li 2 (l/a;) = — — - log 2 x — in log x x > 1. (D-14) 

Sometimes Li 2 (l — x) is referred as Spence function. 

Dilogarithms are a special case of the class of polylogarithms, defined as: 

^HE^S/'*^p. (D15) 

k=l 
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Appendix E 



Harmonic Polylogarithms 



The occurrence of dilogarithms in one loop calculation in QED and QCD is well 
known since a long time. It is also well known that two loops calculations require 
the introduction of the Nielsen's generalization of dilogarithms, also known as poly- 
logarithms. 

Nielsen's polylogarithms are defined as 



It is known that, while the class of Nielsen's polylogarithms is sufficient for a com- 
plete analytical calculation of two loops amplitudes, it is not when the number of 
loops becomes higher or several different scales are involved in the calculations. 
In these case a generalization of Nielsen's polylogarithms which allows a complete 
analytical calculations is provided by the class of harmonic polylogarithms (hpl's) 
introduced in jZHj; by Remiddi and Vermaseren. 

Following the formalism in |75j . the harmonic polylogarithms are labelled by a set of 
arguments contained in the if- dimensional vector m w and are indicated as H(m w ; x). 
For example for w — 1 they read 




(E.l) 



H(0;x) 



log x, 




(E.2) 



Their derivatives are: 



—H(a;x) = f{a;x), 
where a = 0, 1, —1 and the functions f(a; x) are 



(E.3) 



/(0;aO 



1 



x 



/(i;z) 



l 



1 — X 
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/(-l;*) = j^. (E.4) 



One could start by defining the functions in (|E.4|I as the basis to build harmonic 
polylogarithms of weight 1 or higher. 
In fact higher weight hpl's are defined as 

H(m w ;x) = / dx' f(a;x')H(m w -i;x'), (E.5) 
Jo 

where 

m w = (a,rh w ^i), (E.6) 



being a the leftmost index of m w . The definition in (|E.5|) is valid for 7^ 0, that 
is if almost one index is different from zero. 
For rh w = W 

H(d w ,x) = ±- log- x. (E.7) 
W ! 

In compact form the derivative of a general harmonic polylogarithm of weight w is 

-^-H(rh w ; x) = /(a; x)H{m w ^ x). (E.8) 
ax 

For w = 2 the harmonic polylogarithms can be expressed in terms of well known 
functions as logarithms and dilogarithms: 



#(0,0; 


;r) 


= ^log 2 x, 


#(0,1; 


x) 


= Li 2 (x), 


#(0,-1; 


x) 


= -Li 2 (-x), 


#(1,0; 


x) 


= —log xlog(l — x) + Li 2 (x), 


#(1,1; 


x) 


= 2 log2 ( 1-:r )' 


#(i,-i; 


x) 


1 — X 

= Li 2 { ) log21og(l x) 


#(-1,0; 


x) 


= logx log(l — x) + Li 2 (— x), 


#(-i,i; 


x) 


= L% 2 ~ lo § 2 log(l + x) 


#(-i,-i; 


x) 


= ilog 2 (l + x). 



Li, ( \ 



(E.9) 

For polylogarithms of weight w = 3 the same thing happens, that is they can be 
all expressed in terms of logarithms, dilogarithms and Nielsen's polylogarithms of 
weight 3. For higher weights this is no longer true: the class of harmonic polyloga- 
rithms in this case is a wider class with respect to Nielsen's polylogarithms. 
The behaviour at the border of the integration domains is the following: 



if m w = 0, trivially H(0 W ; x) has a logarithmic divergence for x — > (see (jE.7l 
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if rh w 7^ 0, H(rh w ; 0) = for x = 0; 

if m w 7^ 1, H(rh w ; 1) is finite; 

if m w = 1 and m w _i = 6 w -i, H(rh w ; 1) is finite; 

if m w = 1 and m w _i 7^ w -i, H(rh w ; 1) has a logarithmic divergence in x — 1, 
and the term with the highest divergence is log p (l — x), where p is the number 
of the leftmost indices equal to I. 
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Transverse Momentum Distributions in B Decays 
by Roberto Sghedoni 

In the last few years b physics has become a very important topic in parti- 
cle physics: for the first time colliders dedicated to the study of the b quark, 
BELLE and BABAR, have been built and have begun to give precision results. 
The quark beauty has a peculiar role in the Standard Model: at present, this the- 
ory describes the interaction of elementary particles via electroweak and strong 
interactions. For several reasons the Standard Model is not believed to be the 
final theory of elementary particles and efforts have been done to find processes 
confuting its predictions. 

b physics is very sensitive to the limits of the Standard Model: it involves many 
parameters of the theory, poorly known before the construction of the b fac- 
tories, whose precise knowledge can confirm or violate theoretical predictions. 
Important parameters studied are the mass of the b quark and CKM matrix 
elements, for example. 

The rich phenomenology of the physics of the beauty quark is dominated by the 
transition 

b^c 

which is overwhelming in almost the whole phase space. 

Beside the dominant transition, processes with a lower branching ratio are 
present: they are called rare decays and are suppressed for several reasons. 
For example the transition 

b — > u + e~ + v 

is Cabibbo suppressed, because the Cabibbo matrix element involved in the 
process, V u b, is very small compared to V c b- V u b is one important parameter to 
study to constraint the Standard Model and to find possible violations. 
Another important rare process is the so called radiative decay 

b — > s-f 

which is the main topic of my thesis. 

It has been study since the beginning of the 90s for several reasons: 

• it is mediated by a flavour changing neutral current: this can happen in 
the Standard Model only by loop mediated processes and this is the reason 
of its suppression; 

• because of the previous point it is very sensitive to new physics effects, 
such as the presence of a new Higgs doublet or supersimmetric particles, 
though present experimental data does not suggest this eventuality; 

• it involves non perturbative functions which arc universal and can be 
extracted by this process to be used into other processes such as b — > 
u + e~ + V; 
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• it shows a non trivial interplay between strong and electroweak interac- 
tions which permits the application of advanced method of field theory, 
such as effective theories. 

In particular the last point has been widely studied since many years, building 
an effective hamiltonian which describes this process and this approach has 
been applied to calculate the inclusive rate and the photon spectrum by several 
authors. 

We decided to study the transverse momentum distribution of the strange quark 
with respect to the direction fixed by the photon. This kind of distribution offers 
interesting theoretical features: 

• like the photon spectrum it shows the appearance of large logarithms 
which need to be resummcd, according to the well known theory of rcsum- 
mation, applied in QCD in high energy processes involving light quarks. 
In this case we test this approach in a lower energy process involving an 
heavy quark; 

• the resummed distribution shows peculiar singularities which signals the 
appearance of non perturbative effects. As we stated the treatment of 
perturbative effects, which is central in QCD, cannot be performed as 
previously done in literature for the photon spectrum; 

• the fixed order contribution, needed to get a resummed distribution with 
next-to-leading accuracy, has been performed by mean of harmonic poly- 
logarithms, according to a technique introduced by Remiddi and Ver- 
maseren; 

Finally we obtain a resummed distribution with next-to-leading accuracy, which 
can be compared with future data. 

To conclude, we underline that, even if this a rare process, we decided to study 
it because it has a very simple kinematics and it should be a good way to check 
the reliability of the perturbative expansion in b physics. 

Since its discovery the perturbative expansion has been widely applied to b 
physics, in inclusive quantities, such as the inclusive rate, and semi-inclusive 
quantities, such as the photon spectrum, where the improved perturbation the- 
ory, applying the resummation, has become a standard technique. 
However, due to the uncertainties of many parameters, it is not so clear if the 
perturbative expansion is reliable at so low energies (m b ~ 5 GeV) and, a 
fortiori, this matter is valid for semi-inclusive distribution as the one we have 
considered. 
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